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1.  Introduction 


Numerous  sensing  modalities  can  produce  3D  data  sets.  Two  very  important  surveillance  modalities  are  video 
(2D  spaee  plus  time)  and  3D  Inverse  Synthetic  Aperture  RADAR  (1SAR).  It  is  the  goal  of  this  task  to  examine 
these  specific  modalities  and  1)  for  video  -  determine  fundamental  bounds  on  sub-sampling,  resolution  and 
reconstruction  such  that  large  video  data  volumes  ean  be  moved  through  communication  channels;  and  2) 
formulate  and  demonstrate  a  3D  1SAR  image  reconstruction  proeess  to  be  used  for  penetrative  reconstruction  of 
river  boat  internals  as  they  float  by  a  small,  low  frequency  RADAR  The  envisioned  system  eould  provide 
remote,  unattended  river  intelligence,  surveillance  and  reconnaissance  (1SR)  for  several  months  in  a  concealed 
package  located  near  tropical  rivers,  mountain  passes,  and  locations  of  special  interest. 


The  kernel  of  this  projeet  is  the  development  of  reconstruction  algorithms.  Specifically  the  development  and 
implement  of  novel  algorithms  for  generating  3D  representations  of  targets  and  their  contents  using  a  single 
penetrative  RADAR  system.  The  application  area  for  this  technology  is  river  boat  surveillance.  An  overview 
figure  of  this  application  area  is  shown  in  Figure  1. 


Penetrative  3D  Riverboat 
Imaging 


Objectives 


Ovfi  all  Small,  self-sustaining  tree 
based  system  that  can  rapidly  image 
the  internal*  ofnver  boat*  via  their 
natural  motion  down  a  nver  revealing 
weapons,  contraband,  etc  . 

Study  Provide  proof -of -concept  that 
penetrative,  3D,  fine  resolution  ISAR 
reconstruction  can  occur  using  an 
unattended,  low  power,  low 
frequency  Radar  on  full  scale 
Riverboat  traffic 


Figure  1:  The  algorithm  development  is  in  support  of  a  3D  penetrative  imaging  system 

which  monitors  boat  traffic. 


The  concept  of  operations  for  this  sensing  modality  places  restrictions  on  the  type  and  amount  of  diversity  that 
ean  be  used  in  3D  reconstruction.  Because  the  sensor  is  surveilling  a  specific  part  of  a  river  and  is  fixed  in  a 
tree,  the  boat  must  provide  the  angular  diversity  necessary  for  3D  reconstruction.  We  will  not  consider  an 
interferometrie  system  due  to  packaging  constraints.  Figure  2  shows  where  a  potential  system  will  be  placed. 
The  implications  of  these  placements  is  that  the  boat  will  cither  be  moving  straight  through  the  RADAR  beam 
as  it  moves  down  a  straight  segment  of  the  river  or  will  seribe  a  eurve  as  it  move  around  a  river  bend.  We  are 
exploring  the  amount  of  3D  information  inherent  in  that  type  of  motion 
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The  ideal  types  of  motion  include  roll/piteh/and  yaw  which  fills  a  3D  measurement  space  with  data  However, 
this  isn’t  compliant  with  the  eventual  deployment  of  sueh  a  system  and  would  motivate  installing  the  RADAR 
to  perform  surveillanee  only  on  areas  of  the  river  that  have  rapids  which  is  where  boats  don’t  wish  to  go.  The 
straight  and  eurved  trajectories  give  less  3D  fill  Specifically,  as  the  boat  moves  through  the  large  RADAR 
beam,  there  is  a  diversity  in  azimuth  angle  and  in  depression  angle  which  provides  some  inherent  3D  data  fill. 
We  will  be  concentrating  on  that  seenano  in  our  development  and  testing. 


Figure  2:  The  placement  of  the  RADAR  along  a  river  constrains  the  amount  of  angular 
diversity'  that  is  created  due  to  the  boat  motion. 

Because  the  overall  sensing  processing  goal  is  penetrative  imaging  of  river  boats  to  determine  their  internal 
structure  and  whether  they  are  carrying  contraband,  the  sensing  modality  is  a  low  frequency,  high  bandwidth 
RADAR.  We  are  operating  in  the  500MHz  to  >  2GHz  frequency  bands  to  allow  for  penetration  with  modest 
RADAR  powers  along  with  sufficient  bandwidth  to  isolate  specific  scattering  centers  in  range  sueh  that  they 
can  be  traek  over  time.  This  quarterly  report  outlines  the  progress  that  has  been  made  in  the  algorithm 
development  and  also  describes  the  supporting  data  collections  and  simulations  used  in  the  algorithm 
development. 

We  have  performed  signature  simulations  and  data  collections  to  provide  input  to  be  used  by  the  algorithm 
development  Figure  3  shows  a  2D  SAR  example  based  on  a  boat  simulation.  The  1SAR  reconstruction 
algorithm  is  a  phase  estimation/eompensation  algorithm  on  specific  scattering  center  tracks  using  a  rigid  body 
assumption  of  motion  that  ties  the  various  tracks  together  spatially  Invariant  geometrical  measurements, 
resolution  of  three-dimensional  shape  via  comparison  of  multiple  views,  fusion  with  optical,  more  literal 
interpretation,  simple  relation  to  the  user’s  point  of  view,  greater  target  discrimination  capabilities. 

This  report  is  organized  as  follows:  Section  2  provides  an  overview  of  the  processing  approach  and  algorithm 
developments  attained  during  this  reporting  period.  Section  3  describes  the  simulations  performed  and  provides 
example  of  the  simulations  that  were  used  in  the  signature  analysis  and  algorithm  development.  Section  4 
discusses  the  local  collections  that  were  made  with  the  MTR1  RADAR  assets  to  provide  empirical 
measurements  for  analysis  and  algorithm  development  Section  5  discusses  the  data  and  provides  results  of  the 
data  analysis.  Section  6  provides  information  in  regards  to  the  antenna  that  was  designed  for  this  task  and 
diseusses  its  performance,  A  summary  is  then  provided  in  Section  7. 
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Xpatch  Simulation  with  Weapons  Cache 
Radar  image  formed  using  wavenumber  migration 


However  Need  to  Determine 


Weapons  Inside  and  Outside 

Cache  Inside 

Canal  Boat,  With  Weapons  Cache 


Stove  Pipe 


Front  Exterior 
Bulkhead 


Engine-^ 

View  of  model  interior 


-  Bandwidth.  1  7GHz  0  3-2  GHz 
-Azimuth  extent*  *141  deg 
113  degrees  elevation 

Added  2  x  1  x  0.5  meter  metal 
cuboid  in  center  of  hold 


Rear  Exterior 
Bulkhead 


16 


2t 


11 


It 
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Figure  3:  Penetrative  RADAR  and  2D  reconstruction  produces  signatures  that  can  be 
interpreted  as  revealing  boat  infrastructure  and  cargo. 


2.  Overview  of  Processing  Approach 

The  flow  diagram  shown  in  Figure  4  presents  a  very  high  level  concept  for  a  processing  chain  for  the  radar 
signals.  It  has  ten  major  components.  For  each  component  of  the  diagram,  the  envisioned  system  will  contain 
multiple  subcomponents,  and  there  are  many  possible  variations  in  the  implementation  details.  Later  diagrams 
will  illustrate  some  of  those  possible  components.  But  first  we  will  describe  the  purposes  of  the  major 
components. 
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Figure  4  High  level  How  diagram  For  a  possible  processing  system 


The  First  component  comprises  an  initial  data  processing  stage  in  which  many  oF the  specific  radar  and 
environment  dependent  eharaeteristies  oF  the  collected  signals  may  be  addressed  These  issues  include 
conversion  oF  the  data  From  w  hatever  Format  is  most  convenient  For  data  collection  into  a  Format  which  is 
convenient  For  processing.  Other  issues  include  mitigation  oF  any  deterministic  signals  in  the  data  which  have 
been  introduced  by  electrical  leakage  in  the  collection  hardware.  This  appears  to  be  a  pervasive  problem  oF 
small  compact  systems.  This  processing  stage  may  also  need  to  mitigate  undesirable  Features  in  the  impulse 
response  oF  the  collection  system,  to  mitigate  interFerenee  signals  w  hich  were  injected  into  the  data  by  nearby 
sources  oF  radio  Frequency  signals,  mitigate  eFFeets  oF  gaps  in  the  set  oF  collected  radio  Frequencies  due  to 
transmit  restrictions,  and  handle  irregularities  in  the  collection  times  oFthe  radar  sweeps. 

The  seeond  component  oFthe  envisioned  processing  system  achieves  at  least  three  major  Functions.  First,  the 
target  signals  contained  in  the  data  are  modified  by  the  processing  so  that  reductions  are  achieved  in  the 
apparent  range  extent  oF  the  target  signals,  the  maximal  speed  (opening  or  closing  rate)  oFthe  target  signals,  and 
the  accelerations  oF  the  target  signals.  Seeond,  these  reductions  enable  a  compression  oF  the  target  signals  into 
Few  digital  samples,  which  makes  enables  more  eFficient  subsequent  processing  (and/or  signal  exfiltration). 
Third,  the  signal  to  clutter  plus  noise  ratio  is  enhanced  in  the  digital  remaining  samples. 

The  third  component  oF  the  processing  begins  the  automated  analysis  oFthe  target  signals.  Its  purpose  is  to 
separate  signals  From  distinct  parts  oF  the  target,  and  to  estimate  the  changing  ranges  to  eaeh  oF  those  parts  oF 
the  target.  This  stage  is  needed  because,  encoded  within  the  extracted  pulse  by  pulse  descriptions  oF  the  ranges 
to  the  multiple  distinct  parts  oF  the  target,  is  all  the  information  whieh  the  radar  signals  contain,  concerning  the 
size,  shape,  and  relative  motions  (angular  motions)  oF  the  observed  parts  oFthe  target.  This  processing  stage 
produces  range  histories,  For  a  set  oF  distinet  physical  components  oF  the  target.  Eaeh  range  history  is  a  Function 
oF  time,  sampled  at  the  radar’s  pulse  (or  sweep)  rate,  whieh  give  the  distance  From  the  radar  to  a  particular 
component  oF  the  target.  This  stage  oF  the  processing  extracts  the  range  histories  For  as  many  target  components 
as  can  be  accurately  Followed,  but  it  does  not  need  to  extract  range  histories  For  all  oF  the  target  components.  1  F 
the  extracted  range  history  estimates  are  sufficiently  precise.  Four  distinet  histones  will  suffice  For  the 
subsequent  processing  stages,  but  more  will  generally  be  extracted  to  provide  redundancy  whieh  subsequent 
components  may  use  to  mitigate  residual  errors  in  the  estimates.  At  this  stage  of  the  processing  the  system  does 
not  require  any  knowledge  of  which  part  of  the  target  a  particular  range  history  estimate  is  associated  with.  The 
association  of  the  range  histories  with  specific  physical  components  of  the  target  will  be  determined  later,  by 
subsequent  stages  of  the  signal  processing. 

The  Fourth  component  of  the  processing  system  examines  the  range  history  estimates  jointly.  This  is  a 
robustifieation  stage.  This  stage  is  not  needed  if  all  of  the  range  estimates  From  the  previous  stage  are  good 
estimates  oFthe  ranges  to  distinet  components  of  the  target.  But  we  include  it  the  system  design,  because  we 
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have  developed  methods  to  use  geometric  constraints  on  the  joint  behavior  of  the  true  range  histories  to  rank  the 
quality  of  a  set  of  estimated  range  histories,  and  to  eliminate  estimates  which  may  contain  large  errors.  (There 
arc  various  possible  ways  in  which  an  erroneous  range  history  may  have  been  obtained.  Imagine,  for  example, 
that  the  system  has  extracted  the  range  history  from  a  reflection  from  a  log,  which  happened  to  float  through  the 
scene  while  it  was  being  traversed  by  the  target.) 

The  fifth  component  uses  the  range  histories  to  extract  the  three-dimensional  si/c  and  shape  of  the  configuration 
of  the  target's  reflective  components.  It  extracts  a  complete  three-dimensional  description  of  the  geometry  of  the 
configuration  of  the  set  of  components  that  reflected  the  radar  signals  (from  w  hich  the  range  histories  were 
extracted).  All  characteristics,  such  as  distances,  angles,  areas,  and  volumes,  of  the  configuration  become 
known,  and  each  of  the  range  histories  becomes  associated  with  a  specific  location  on  the  target  The  processing 
strategy  for  this  part  is  based  on  geometric  theorems  which  describe  how  the  measured  ranges  to  distinct  points 
on  the  target  are  constrained  by  the  size  and  shape  of  the  target.  These  constraints  lead  to  systems  of  equations 
which  determine  a  fundamental  set  of  geometric  invariants,  w  hich  in  turn  determine  the  size  and  shape.  At  this 
stage  of  the  processing,  some  of  the  measurements  of  target  and  cargo  characteristics,  which  might  be  directly 
utilized  by  an  intelligence  analyst,  begin  to  emerge  from  the  processing.  Yet,  there  is  more  processing  to  be 
done,  so  as  to  extract  additional  information  from  the  data,  and  to  put  that  information  into  useful  forms  for  the 
end  user. 

The  sixth  component  uses  the  detailed  size  and  shape  information,  from  the  previous  processing  stage,  along 
with  the  associated  range  histories,  to  determine  all  of  the  observable  motions  of  the  target.  Because  the 
previous  stage  has  determined  the  shape  of  the  configuration  of  the  reflectors,  and  the  range  histories  keep  track 
of  the  pulse  by  pulse  distances  from  the  reflectors  to  the  radar,  constructing  the  pulse  by  pulse  motions  becomes 
just  a  matter  of  solving  a  small  set  of  equations  at  each  pulse  time.  The  motion  solution  yields  the  specific 
direction  from  which  the  target  was  illuminated  at  each  pulse  time.  It  also  refines  the  initial  range  to  target 
estimates,  and  enables  the  calculation  of  the  pulse  by  pulse  range  to  any  point  (reflective  or  not)  on  the  target. 
The  seventh  component  uses  the  motion  information  to  form  various  image  products.  These  products  arc  based 
on  all  the  information  provided  by  the  sensor  signals.  The  motion  estimates  from  the  previous  stage,  provide  all 
of  the  auxiliary  information  which  is  needed  by  conventional  SAR/1SAR  image  formation  processing 
algorithms  The  image  formation  processing  applies  the  motion  dependant  transformations  to  the  clutter  reduced 
radar  signals,  which  were  provided  by  the  second  stage  of  the  processing.  Conventional  two-dimensional 
images  can  be  formed  from  subsets  of  the  radar  signals,  using  conventional  methods  (polar  formatting, 
migration,  chirp  sealing,  two-dimensional  back-projection, ...)  Overlapping  subsets  can  be  used  to  produce 
movies  ("dynamic  images")  showing  the  motions  of  the  target  (and  possible  indications  of  motions  of  objects  or 
persons  on  the  target)  Also  all  the  collected  pulses  can  be  simultaneously,  to  produce  images  with  three  spatial 
dimensions  Conventionally,  this  is  done  with  baek-projeetion,  but  more  advanced  methods  can  also  be  used  to 
produce  better  images  by  exploiting  the  natural  sparseness  of  illuminated  reflective  surfaces  in  three- 
dimensional  space. 

The  remaining  three  components  arc  not  yet  as  well  defined  by  our  current  work,  but  arc  included  in  the 
description  because  of  their  importance  for  expanding  the  audience  for  the  products  of  the  system 

The  eighth  component  is  planned  to  enhance  the  utility  of  the  shape  measurement  and  image  products  by 
corresponding  the  products  with  stored  additional  sources  of  target  and  cargo  knowledge,  such  as  models  and 
images  of  known  boat  types,  and  previously  collected  data  sets. 

The  ninth  component  will  use  knowledge  of  boat  design,  cargo  characteristics,  boat  bcha\ior  patterns  and 
regional  knowledge,  together  with  the  model  and  image  correspondences  from  the  previous  component,  to 
enhance  the  interpretability  of  the  sensor  products. 

The  tenth  envisioned  component  is  conceived  as  primarily  an  information  formatting  and  display  component. 
Having  done  everything  that  can  be  done  efficiently  by  an  automated  system  to  prepare,  the  last  component 
must  be  designed  to  facilitate  the  exploitation  of  the  combined  information  by  the  user. 
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3.  Simulations 


Various  data  simulations  were  eonstrueted  during  the  course  of  the  program.  T  hese  simulations  supported 
algorithm  development  both  by  providing  test  data  in  the  period  before  actual  collected  data  were  available,  and 
to  simplify  and  separate  phenomenology  for  closer  study. 

Two  types  of  simulations  were  performed,  one  consisting  of  simple  point  seatterer  responses,  the  other  utilizing 
the  Xpatch  simulation  tool.  The  simple  point  seatterer  simulations  run  quickly,  providing  simulations  quickly. 
The  Xpatch  simulations  give  some  insight  into  the  interactions  due  to  penetration  through  various  dielectric 
materials,  as  well  as  multi-bounce  effects.  Both  are  valuable  tools  during  algorithm  design  and  development. 
To  support  the  simulations,  we  obtained  a  faceted  model  of  a  river  boat.  We  assigned  materials  to  the  boat  to 
match  the  project  eon-ops:  mainly  wood  construction,  with  various  smaller  metal  parts.  The  boat  was  also 
given  a  motor.  Figure  5  shows  the  canal  boat  model  used  in  the  simulations. 

Metal  smokestack, 
rail,  and  motor 


_  and  hull 

Figure  5  Canal  Boat  Model 

T  he  boat  was  placed  in  a  plane  to  simulate  being  on  the  surface  of  water. 


Figure  6  Canal  Boat  with  Hater  Plane 
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The  v  arious  simulations  were  designed  to  emulate  the  parameters  of  the  expected  collected  datasets. 

Simulations  were  done  in  the  near  field  at  low  frequencies  (300-2000  MHz)  to  simulate  penetration  of  RADAR 
signals  through  the  dielectric  materials. 

Canal  Boat  7  Point  Scattcrcr  Simulations 

To  help  us  understand  issues  related  to  both  near  field  effects  and  scattcrcr  interaction  during  typical  boat 
motion,  we  created  a  7  point  model  based  on  the  canal  boat.  We  chose  the  7  points  to  span  the  3  dimensions  of 
the  boat,  to  support  simulation  of  typical  scattcrcr  responses  and  interactions,  and  provide  an  unambiguous 
signature.  Figure  7  shows  the  locations  of  the  7  points  on  the  model  (note:  The  model  is  being  shown  partially 
transparent  to  illustrate  scatterer  locations). 


Figure  7  Canal  Boat  7  Point  Model  Locations 

The  points  in  the  7  point  model  were  given  the  motions  expected  of  a  canal  boat  going  around  a  bend  in  a  river. 
Figure  8  shows  some  points  around  a  river  bend,  and  a  possible  sensor  location,  that  were  used  to  generate  a 
path  for  the  points. 
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Figure  8  River  Bend  Simulation  Motion 


This  produced  the  types  of  ranges  histories  we  expected  to  see  in  the  collected  data  Figure  9  shows  the  ranges 
to  the  points  in  the  7  point  model.  Figure  10  shows  the  range  compressed  signals  created  from  the  simulated 
signal  history  created  from  those  ranges,  and  F  igure  1 1  show  s  a  closeup  of  the  range  compressed  signals  during 
the  period  when  the  boat  crosses  broadside  to  the  RADAR,  causing  the  scatterers  to  cross 
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Ranges  to  Points  in  7  Point  Model 
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Figure  9  Ranges  to  Points  in  1  Point  Canal  Boat  Model 


Range  Histories  From  Simulated  signal  History  of  7  Point  Model 
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Figure  1 0  Simulated  Range  compressed  signals  of  7  Point  Canal  Boat  Model 


Three  Dimensional  Imaging  Studies  and  Demonstration 


MTRI  •  1-9 


Sweep  Number 

Figure  J I  Closeup  of  Simulated  Range  compressed  signals  of  7  Point  Canal  Boat  Model 


Xpatch  Canal  Boat  Simulation 

An  Xpatch  simulations  was  run  to  provide  insight  into  the  effects  of  dielectric  materials  and  multi-bounee 
effeets  on  the  range  compressed  signals.  The  eanal  boat  with  water  plane  was  used  as  the  target,  and  it  was 
given  the  same  motion  as  the  7  point  model.  We  also  used  the  same  RADAR  parameter.  The  range 
compressed  signals  for  the  eanal  boat  Xpatch  model  arc  shown  in  Figure  12  and  Figure  13. 


Sweep  Number 

Figure  12  Simulated  Range  compressed  signals  of  Xpatch  Canal  Boat  Model 
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Figure  13  Closeup  of  Simulated  Range  compressed  signals  ofXpatch  Canal  Boat  Model 


Foam  Cart  Point  6  Point  Scatterer  Simulation 

In  preparation  for  our  data  collection,  wc  performed  some  local  collections,  which  arc  discussed  later  in  this 
document.  One  of  these  collections  was  a  cart  covered  in  foam,  w  ith  6  reflectors  mounted  on  the  top.  To  help 
us  understand  the  phenomenology,  we  used  the  point-to-point  distances  between  seatterers  to  create  a 
simulation  of  this  cart,  and  gave  it  a  straight  line  motion  in  front  of  the  sensor.  The  locations  of  the  seatterers 
arc  shown  in  Figure  14 
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Figure  14  Point  Seatterers  for  6  Point  Foam  Cart  Model 
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We  gave  the  points  the  motion  shown  in  Figure  15,  and  these  produced  the  range  histories  shown  in  Figure  15. 
In  Figure  16  and  Figure  17  the  range  compressed  signals  produced  from  the  simulated  signal  history  arc  shown. 
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Range  Histories  From  Simulated  Signal  History  of  6  Point  Model 
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Figure  16  Simulated  Range  compressed  signals  of  6  Point  Foam  Cart  Model 
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Ranges  to  Points  in  6>  Point  Foam  Cart  Model 
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Range  Histories  From  Simulated  Signal  History  of  6  Point  Model 
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Figure  1 7  Close  up  of  Simulated  Range  compressed  signals  of  6  Point  Foam  Cart  Model 


4.  Local  Collections 

Prior  to  field  collections,  it  was  prudent  to  characterize  the  proposed  system  configuration  under 
controlled  conditions,  and  in  an  environment  in  which  rapid  changes  and  corrections  could  be  made 
as  needed.  To  this  end,  we  mounted  the  RADAR  system  on  the  roof  of  our  building  with  a  view  of  a 
clear,  grass-covered  area  that  had  a  geometry  similar  to  proposed  field-collections. 

The  system  itself  consisted  of  an  Akela  radar,  the  McCorkle  antennas,  a  1 W  power  amplifier,  and  20ft 
low-loss  cables.  The  radar  specification  are  depicted  in  Figure  18  below: 
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Akeia  Radar  Specifications 


380MHz  -  3000MHz  FMCW 
Sample  rates  from  3KHz  -  90KHz 
Bui  ft -in,  selectable,  low-nos#  amplifier 
Built-in  range-gating 
IX  putec  length  24m  -  1556ns 
RX  delay  18m  1556ns 
RX  gala  18n>  -  I556na 
Interpulse  delay  1 8ns  1 566ns 
Output  power 

1  /dBm  typical  20d8m  mexirmm 
Input  signal  limits 

♦  13dBm  (LNAoff) 

-7  dBm  (LNAon) 

4  -  chann#l 
Operates  -n  RX/TX  pairs 

2  position  encoder  inputs 


Figure  1H:  Specifications  of  the  Akeia  RADAR  system  used  for 

collections. 

The  geometry  of  the  collection  is  depicted  in  Figure  19  below: 


Measurement  Setup 
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Figure  7  9:  Local  collection  geometiy. 
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The  overhead  view  above  does  not  accurately  reflect  the  vehicles  that  were  present  during  the 
collection.  The  range  gate  of  the  radar  was  adjusted  to  reject  the  vehicles  present. 

Targets  were  then  built  and  mounted  on  a  large  piece  of  construction  foam.  The  foam  itself  is 
transparent  at  the  frequencies  of  interest,  but  served  as  a  rigid  platform  on  which  the  targets  could  be 
arranged  as  needed.  The  targets  consisted  of  trihedrals,  spheres,  top-hats,  and  an  RPG  in  its  wooden 
shipping  crate.  Some  target  arrangements  are  shown  in  Figure  20  below.  The  RPG  used  in  the 
collection  is  shown  in  Figure  21. 


F«r  Range 
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Figure  20:  Target  layout  used  on 


Photos 


FuN  Cart  Configuration  Equipment  Setup 


cart  for  local  collections. 


The  cart  itself  had  30m  ropes  attached  at  the  comers  to  allow  it  to  be  dragged  by  staff  that  were 
outside  the  range-gate  extent.  The  cart  w  as  started  to  the  right  of  the  radar  beam,  and  dragged  at  a 
relatively  constant  rate  (estimated  to  be  0.75m/s)  through  the  beam,  and  well  beyond,  to  provide  a 
target-free  background  estimate  in  every  data  set  collected. 


Figure  21:  RPG  used  as  target.  RPG  was  in  a  wooden  shipping  crate  during  collections. 
Experiments 


Data  were  collected  under  these  conditions: 

1.  5  00 M  H  z-3  000M  H  z 

2.  501  frequency  samples  evenly-spaced  in  frequency  over  this  range 

3.  700mW  estimated  effective  radiated  power 

4.  Various  target  arrangements,  ranging  from  an  empty  cart  (with  mounting  hardware, 
shown  in  the  figure  above),  to  a  crowded  arrangement  of  up  to  6  targets.  The 
experiments  performed  were: 


Three  Dimensional  Imaging  Studies  and  Demonstration 


MTRI  •  1-15 


1 .  Empty  cart 

2.  1  trihedral 

3.  2  trihedrals 

4.  all  targets  (placed  as  indicated  in  the  figure  above,  left) 

5.  RPG  on  foam,  angled  away  from  radar 

6.  RPG  on  foam,  angled  toward  the  radar 

7.  RPG  elevated  above  foam,  angle  away  from  radar 

8.  RPG  elevated  above  foam,  angled  toward  the  radar 

9.  RPG  in  wooden  box  under  the  variations  indicated  above 

The  reason  for  varying  the  RPG  elevation  was  to  be  able  to  assess  the  effects  of  having  targets  very 
close  to  a  ground-plane.  It  was  felt  that  this  might  mimic  targets  being  in  a  boat  on  the  bottom  of  the 
hull  (close  to  the  water),  as  well  as  a  target  on  a  deck  above  the  water. 

Example  Results 

Shown  below  are  some  of  the  quick  look  2D  results  obtained  during  these  tests.  T  he  data  analysis 
section  will  describe  the  3D  processing  of  these  data. 

The  first  example,  shown  in  Figure  22,  is  a  range-compressed  image  of  the  foam  cart  with  2  trihedrals 
present.  All  the  range-compressed  images  shown  are  formed  by  taking  the  Fourier  Transform  of  each 
frequency  sweep,  and  displaying  the  time  course  of  the  sweeps.  Increasing  time  is  to  the  right,  and 
increasing  range  from  the  radar  is  towards  the  top  of  these  images. 

R*ng«  Profile  from  roof  r*of_1  (OIFFED,  T©p40d9) 


Figure  22:  Range  compressed  signature  of  two 
trihedral  configuration  as  the  sled  moves  by  in  time. 

As  the  targets  enter  the  radar  beam,  they  become  increasingly  bright,  their  ranges  decrease  until  they 
are  broadside  to  the  radar,  and  then  they  move  through,  and  away  from  the  beam.  This  image  was 
processed  with  background  clutter  suppression,  and  the  target-to-clutter  ratio  is  seen  to  be  over  30dB. 
This  would  be  considered  ideal  data  for  the  analysis  we  wished  to  perform. 
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Another  example  showing  the  range -compressed  image  for  the  sled  configuration  using  all  targets 
present  as  was  shown  in  Figure  20  is  shown  in  Figure  23. 


Range  Profck  Irem  roof  rx>Ml_l2J0_»N>rget2_data  (OF FED,  Top40dB) 
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Figure  23:  Range  profile  of  sled  with  full  target 
complement. 

Again,  these  data  appear  to  be  of  high  quality;  there  are  few  system  artifacts  present  after  clutter 
suppression,  and  the  dynamic  range  of  the  targets  is  over  30dB. 

Although  we  did  not  have  motion  information  recorded  for  these  data,  it  was  observed  that  the  staff 


Two  Targets 


Figure  24:  Example  2D  ISAR  image  of  two  trihedral 
configuration. 
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were  able  to  drag  the  cart  at  a  very  regular  pace,  allowing  fully  2-D  ISAR  images  to  be  formed 
without  the  need  for  adaptive  phase  estimation.  Because  the  target  was  in  the  near-ficld  of  the  radar, 
and  wavefront  curvature  was  significant,  we  used  a  range-migration  algorithm  to  form  images. 

An  example  of  such  an  image  is  show  n  in  Figure  24  for  the  case  of  2  targets  mounted  on  the  foam: 

In  these  images,  increasing  range  from  the  radar  is  to  the  right,  but  due  to  a  sign  error  in  the 
processing  code,  the  image  is  the  mirror  of  the  actual  geometry.  However,  it  can  be  seen  that  the 
target  impulse  responses  are  quite  compact  and  symmetric,  indicating  that  the  uniform  linear  path  that 
was  estimated  is  a  good  approximation  to  the  actual  motion.  There  is  at  least  40dB  of  dynamic  range 
in  the  data,  and  systematic  effects  are  minimal. 

Figure  25  shows  the  ISAR  image  of  the  RPG: 


RPG,  in  box 


fpof  11 


1  M  0  M  1  IS  I  n  *  IS: 


Figure  25:  ISAR  image  of  the  RPG  showing  its 
flange  and  front. 

Again,  this  is  the  mirror  image  of  the  actual  arrangement,  but  the  ends  of  the  RPG  are  visible  in  the 
image,  as  is  part  of  the  body,  even  through  the  wooden  box,  as  would  be  expected  at  these 
frequencies. 

Conclusions 

The  system  as  configured  performs  well.  Calculation  show  that  there  is  sufficient  power  to  provide  a 
good  signal  level  at  the  receive  amplifier  for  targets  of  this  RCS  at  the  ranges  expected  in  the  field 
collection.  The  antennas  and  cables  appear  to  be  correctly  matched  over  the  band  of  interest,  and  the 
range-gate  settings  can  be  adjusted  to  receive  targets  in  the  ranges  of  interest,  and  reject  targets 
outside  that  region.  The  system  phase  is  stable,  systematic  errors  are  stationary  over  the  times  that  we 


Three  Dimensional  Imaging  Studies  and  Demonstration 


MTRI.  1-18 


are  likely  to  be  collecting,  indicating  that  clutter  suppression  via  background  subtraction  is  a  viable 
approach.  The  RMA  code  has  been  demonstrated  to  be  substantially  correct,  and  is  being 
incorporated  into  a  “quick  look”  capability,  allowing  rapid  assessment  of  the  data  quality  during  field 
collections.  The  data  are  of  a  high  enough  quality  that  they  can  be  used  in  the  testing  and 
development  of  the  3-D  processing  code. 


5.  Data  Analysis 


In  this  section  we  provide  some  example  processing  results  for  the  simulated  and  locally  collected  data  sets 
described  previously.  These  results  show  the  ability  of  the  processing  algorithm  to  determine  and  track 
scattering  mechanisms  in  3D  such  that  a  3D  spatial  reconstruction  can  occur. 

7  Point  Canal  Boat 


The  7  point  canal  boat  simulation  dataset  was  used  to  study  performance  of  existing  signal  separation 
techniques  with  these  datasets.  The  data  were  first  motion  compensated  (Figure  26)  using  the  range 
to  the  centroid  of  the  7  points  at  each  sweep.  The  goal  of  the  compensation  is  to  slow  the  range  rate  of 
the  signals.  By  compensating  the  signals  from  a  moving  target  to  a  target  fixed  point,  the  translational 
motions  of  the  point  are  removed,  and  the  only  remaining  motions  are  due  to  rotations  of  the  target 
about  the  compensation  point.  This  step  reduces  the  number  samples  needed  to  contain  the  responses 
from  the  target.  There  are  two  reductions  enabled.  One  reduction  is  in  the  number  of  radio  frequency 
samples  that  need  to  be  preserved.  This  is  possible  because  the  range  compression  of  the  signal 
occupies  fewer  range  bins  after  the  motion  compensation.  The  other  reduction  is  in  the  number  of 
times  samples.  This  reduction  is  possible  because  the  Doppler  compression  of  the  target  signals 
occupy  fewer  Doppler  bins  after  the  compensation. 
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Figure  26  Range  Compressed  signals  of  Motion  Compensated  7  Point  Model 
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The  boat  signatures  were  then  extracted  from  the  2D  spatial  domain,  to  help  reduced  the  effects  of 
clutter  and  the  computational  burden.  Figure  27  shows  the  result  of  this  extraction. 


Range  Doppler  Map  of  Reduced  7  Point  Canal  Boat  Model 
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Figure  27  Range  Doppler  Map  of  Reduced  7  Point  Model 

The  goal  of  this  earlier  stage  of  analysis  was  to  produce  range  history  paths  for  the  geometry  analysis 
stage.  To  do  a  good  job  of  extracting  these  range  paths,  the  path  extraction  was  done  in  a  3 
dimensional  cuboid  with  dimensions  time,  range,  and  velocity.  The  2  dimensional  time  frequency 
dataset  used  to  produce  the  range  Doppler  map  shown  above  was  converted  into  a  3  dimensional 
dataset  of  range  velocity  slices  by  taking  overlapping  subsets  of  the  frequency  sweeps.  These  subsets 
of  sweeps  were  converted  to  frequency  Doppler,  then  interpolated  along  Doppler  to  produce 
frequency  velocity.  This  interpolation  produces  lines  of  constant  velocity  in  the  output.  Figure  28 
shows  the  frequency  Doppler  and  frequency  velocity  maps  f  or  the  reduced  7  point  model. 
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Frequency/Velocity  Map  of  Motion  Compensated  7  Point 
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Figure  28  Frequency  vs  Doppler  and  Frequency  vs  Velocity  Maps  of  Reduced  7  Point  Model 


Shown  in  Figure  29  are  three  visualizations  of  the  time-range-velocity  cuboid.  This  is  the  cuboid  that 
will  be  used  in  the  path  finding  process.  These  views  clearly  show  the  paths  of  the  scatterers,  and 
illustrate  that  when  paths  converge  in  range,  they  diverge  in  velocity,  providing  the  leverage  needed 
to  separate  the  scatterers. 
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Figure  29  CTRV  Time,  Range,  And  Velocity  Views 
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6  Point  Foam  Cart 


As  described  earlier,  we  performed  a  local  collection  with  a  cart  covered  in  foam  carrying  6  reflectors 
mounted  on  the  top.  To  help  us  understand  the  phenomenology,  we  used  the  point-to-point  distances 
between  scattering  centers  to  create  a  simulation  of  this  cart,  and  moved  it  in  a  straight  line  motion  in 
front  of  the  sensor.  Figure  30  displays  the  range  compressed  signals  associated  with  the  6  point  foam 
cart  model. 


<D 

-Q 


E 


Q. 

<U 

<U 


5 

IS) 


Range  Compressed  Signals 
From  Simulated  Signal  History  of  6  Point  Model 
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Figure  30  Range  Compressed  signals  of  6  Point  Foam  Cart  Model 

Motion  compensation  was  performed  to  allow  us  to  crop  the  dataset  to  help  reduced  the  effects  of 
clutter  and  lessen  the  computational  burden.  The  motion  compensation  was  performed  using  a  shear 
averaging  algorithm.  The  range  compressed  signals  formed  from  the  motion  compensated  signal 
history  are  shown  in  Figure  3 1 ,  while  Figure  32  shows  both  the  range  compressed  signals  and  range 
Doppler  map  resulting  from  the  spatial  cropping. 
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Motion  Compensated  Range  Compressed  Signals 
From  Simulated  Signal  History  of  6  Point  Model 
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Figure  31  Range  Compressed  Signals  from  Motion  Compensated  6  Point  Foam  Cart  Model 
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Figure  32  Range  Compressed  Signals  and  Range  Doppler  Map ,  Reduced  6  Point  Foam  Cart 
Model 


Using  this  reduced  signal  history,  the  time-range-velocity  cuboid  was  created,  as  shown  in  Figure  33. 
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Figure  33  CTRV Time,  Range,  And  Velocity  Views  of  6  Point  Foam  Cart  Model 


Following  construction  of  the  time -range-velocity  cuboid,  the  paths  of  the  scatterers  were  extracted 
This  extraction  was  done  in  a  way  that  respects  physical  constraints  imposed  upon  changes  in  range 
by  the  current  velocity.  We  also  disallow'  any  velocities  that  are  outside  velocities  of  interest  for  the 
problem  at  hand  These  constraints  do  not  allow'  paths  to  make  transitions  that  are  non-physical. 

The  paths  extracted  from  the  time-range-velocity  cuboid  for  the  6  point  foam  cart  model  are  shown  in 
Figure  34  as  black  lines  overlayed  on  top  of  the  range  compressed  signals. 
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Figure  34  Paths  Through  Time-Range-Velocity  Cuboid  for  6  Point  Foam  Cart  Model , 
Overlayed  on  Range  Compressed  Signals 


The  maximal  difference  between  the  estimated  range  histories  and  the  true  range  histories  was  1 .52 
centimeters  and  the  root  mean  square  distance  was  7.52  millimeters.  Thus  the  differences  are  difficult 
to  see  in  Figure  35  which  shows  the  true  and  estimated  values  together. 
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Figure  35  True  and  estimated  range  histories 

Differences  can  be  seen  more  easily  in  Figure  36,  which  shows  the  true  and  estimated  velocities  of 
the  synthetic  reflectors. 
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Geometric  analysis  was  applied  to  both  the  true  and  the  estimated  range  histories.  The  purpose  of  the 
geometric  analysis  is  to  reconstruct  the  observ  able  components  of  the  shape  and  motion  of  the 
scattering  center  configuration.  Four  versions  of  the  analysis  were  examined.  These  are  the  three  and 
two  dimensional  versions  of  the  near  field,  and  the  three  and  two  dimensional  versions  of  the  far  field 
theory.  The  three-dimensional  near  field  theory  is  the  most  general  theory.  It  makes  no  restriction  on 
the  motion  of  the  target,  and  correctly  accounts  for  spherical  spreading  of  the  radio  waves.  But, 
despite  the  fact  that  it  derives  from  the  most  correct  model  for  the  typical  physical  situation,  current 
algorithms  based  on  this  theory  can  perform  poorly  in  certain  constrained  motion  situations.  These 
problems  arise  when  the  illumination  directions  are  not  diverse  enough  to  keep  the  input  data  well 
removed  from  the  singularity  conditions  associated  with  certain  restricted  types  of  motions. 

The  errors  in  the  estimated  ranges  are  expected  to  introduce  some  error  in  the  geometric 
reconstructions.  Additional  errors  may  be  expected  in  the  estimates  from  the  far  field  estimation 
process  because  it  ignores  wave-front  curvature.  Further,  because  the  true  underlying  motion  was 
linear,  the  three-dimensional  reconstruction  problem  is  singular  for  this  case.  This  should  make  the 
three-dimensional  solution  extremely  sensitive  to  noise  in  the  data. 

The  straight  boat  path  is  an  example  of  a  severely  degenerate  motion  case,  and  this  readily  degrades 
the  estimates  formed  with  the  near  field  algorithm.  This  is  illustrated  by  the  estimates  of  the  scattering 
center  configuration  which  were  obtained  from  the  near  field  algorithm.  These  estimates  are  shown  in 
Figure  37. 
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Figure  37:  True  scattering  locations  (blue)  and  estimated  2D(green)  and  3D  (red)  locations 

based  on  the  invariant  algorithm. 

For  both  the  two-dimensional  (green)  and  the  three-dimensional  (red)  near  field  estimates,  the 
coordinates  are  not  near  to  those  of  the  true  coordinates  shown  in  blue.  The  estimated  shape  is  too 
narrow,  too  close  to  being  collinear,  to  match  the  shape  of  the  true  configuration.  The  poor  estimation 
may  be  attributable  to  the  severity  of  the  over-parameterization  of  the  near  field  invariant  equations  in 
the  degenerate  cases. 

The  far  field  estimates  are  much  better.  A  view  from  a  three-dimensional  display,  shown  in,  Figure 
38,  exhibits  the  true  configuration  (in  blue),  the  three-dimensional  far  field  estimate  (in  red),  and  the 
two-dimensional  far  field  estimate  (in  green). 
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Figure  38  True  and  estimated  reflector  configuration  from  far  field  geometry  analysis.  Truth 
in  blue ,  2d  estimate  in  green ,  3d  estimat  in  red. 

This  particular  view  is  nearly  normal  to  the  plane  of  the  two-dimensional  estimates,  so  the  planarity 
of  the  two-dimensional  estimates  (shown  in  green)  is  not  apparent  (as  it  would  be  from  other  points  of 
view).  In  this  view,  the  estimates  appear  to  match  the  truth  reasonably  well,  despite  all  the  reasons 
why  they  might  not.  The  most  obvious  visible  deviation  is  the  significant  separation  of  the  three- 
dimensional  estimate  (shown  in  red)  near  the  right  side  of  the  image,  from  the  true  location  which 
appears  near  the  far  right  edge  of  the  image. 

The  estimation  of  the  coordinates  enables  the  estimation  of  the  illumination  directions,  where  were 
controlled  by  the  motion  of  the  target.  These  illumination  direction  estimates  can  then  be  used  to 
control  image  formation  processes,  which  convert  the  original  signal  history  into  an  image.  Given  a 
three-dimensional  motion  solution,  both  three-dimensional  and  two-dimensional  image  formation 
process  can  be  used.  The  qualities  of  the  image  products  depends  not  only  upon  the  image  formation 
method,  and  the  accuracy  of  the  motion  estimates,  but  also  on  the  geometric  properties  of  the  set  of 
illumination  directions.  For  good  two-dimensional  images,  finding  large  subsets  of  the  illumination 
directions  which  are  nearly  coplanar  are  best.  For  good  three-dimensional  images,  one  should  use  all 
the  data,  not  a  subset,  and  hope  for  illumination  directions  that  densely  sample  a  large  solid  angle. 

For  this  example  the  estimated  illumination  directions  are  shown  in  Figure  39.  The  blue  curve  shows 
the  angle  of  rotation  about  the  principle  axis  of  rotation  as  a  function  of  the  down-sampled  pulse 
numbers.  The  green  curve  gives  an  indication  of  the  wobbling  of  the  principle  axis  of  rotation. 
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Figure  39  Angular  motion  estimates 

From  the  blue  curve,  we  see  that  there  is  more  than  80  degrees  of  rotation  about  the  principle  axis  of 
rotation.  This  suggest  that  the  data  should  support  fine  resolution  imaging  in  the  plane  that  is 
orthogonal  to  the  principle  axis  of  rotation.  The  green  curv  e  shows  that  there  is  also  some  wobble  of 
the  axis  of  rotation.  Hence  the  estimate  indicates  that  there  is  some  three-dimensional  information  in 
the  data.  If  the  green  curv  e  were  constantly  zero,  the  direction  of  the  estimated  axis  of  rotation  would 
be  fixed.  In  that  situation  the  estimate  would  imply  that  a  two-dimensional  image  could  be  formed 
using  all  of  the  data,  and  that  the  resulting  two-dimensional  image  would  contain  no  height 
dependent  defocusing.  In  that  zero-wobble  case,  there  would  be  no  information  in  the  data, 
concerning  variations  in  the  direction  of  the  rotation  axis.  Because  the  variation  in  the  green  curv  e  is 
not  nearly  as  large  as  the  blue  curve,  we  expect  less  resolution  in  the  direction  of  the  principle  axis  of 
rotation. 


Using  the  three-dimensional  motion  estimates,  we  formed  a  three-dimensional  back-projection 
images  from  the  radar  data.  One  view  of  this  three-dimensional  image  is  shown  in  the  Figure  40.  The 
image  also  show  markers  for  the  estimated  three-dimensional  locations  of  the  reflectors  which  were 
track  though  the  signal  data.  These  markers  were  inserted  after  the  three-dimensional  back-projection 
image  was  formed. 
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Figure  40  Thee  dimensional  hack-projection  image  with  estimated  scattering  center 
locations 


As  expected  there  are  two  directions  of  fine  resolution,  due  to  the  bandwidth  of  the  radar  data,  and  the 
large  angle  of  rotation  about  the  principle  axis  of  rotation.  But  in  the  a  direction  of  the  rotation  axis, 
there  is  poor  resolution  due  to  the  limited  angular  diversity  because  there  is  only  modest  variation  in 
the  estimated  axis  of  the  rotation  during  this  collection. 


Because  of  the  slow  variation  in  the  axis  of  rotation,  reasonably  well  focused  two  dimensional  images 
can  also  be  made  from  subsets  of  the  data.  Such  and  image  is  shown  in  Figure  41 . 
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Figure  41  Polar  formatted  two-dimensional  image  formed  by  using  motion  estimates  from  far 
field  invariant  theoiy 
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Foam  Cart  Data  From  Akela  RADAR 


As  described  earlier,  we  performed  a  local  collection  with  a  cart  covered  in  foam  carry  ing  6  reflectors 
mounted  on  the  top.  This  section  describes  the  processing  of  that  dataset. 


To  help  reduce  internal  system  noise  and  stationary'  clutter  elements,  a  background  subtraction  must 
be  performed  on  the  Akela  data.  During  this  step,  a  principal  components  analysis  was  performed  on 
the  signal  history  data  to  determine  which  parts  of  the  signal  are  stationary.  These  stationary 
components  are  then  subtracted  away.  The  results  of  this  process  on  the  Akela  foam  cart  data  are 
shown  in  Figure  42. 
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Figure  42  Akela  Foam  Cart  Data  Before  and  After  Background  Subtraction 


For  the  purposes  of  range  history  path  extraction,  we  are  only  interested  in  those  sweeps  that  contains 
signals  from  the  target  of  interest.  Those  sweeps  that  contain  no  target  signal  were  eliminated.  The 
result  of  this  extraction  is  shown  in  Figure  43. 
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Akela  Data  Foam  Cart  Range  Compressed  Signals 
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Figure  43  Akela  Foam  Cart  Data  Before  and  After  Target  Extraction 


Motion  compensation  was  then  performed  using  a  shear  averaging  technique  to  slow  down  the  signals  and 
allow  spatial  cropping.  The  motion  compensated  range  compressed  signals  and  range  compressed  signals  and 
range  Doppler  map  of  the  cropped  dataset  arc  show  in 
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Figure  44  Akela  Foam  Cart  Data ,  Motion  Compenstation  and  Spatial  Reduction 
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Due  to  antenna  patterns  imposed  by  the  antenna  used  in  the  foam  cart  collection,  it  was  necessary  to  adjust  the 
amplitudes  of  the  values  in  the  signal  history  domain  to  correct  a  frequency  dependent  falloff.  Additionally,  to 
achieve  higher  resolution,  we  also  used  extrapolation  techniques  to  extend  the  frequency  space,  instead  of 
simple  zero  padding.  These  tw  o  techniques  were  applied  as  part  of  the  conversion  from  the  2  dimensional 
sweep-frequency  data  to  the  3  dimensional  time-frequency-vclocity  data.  Figure  45  shows  the  time-range- 
vclocity  data  before  extrapolation  and  whitening,  while  Figure  46  shows  the  data  after. 
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Figure  45  Akela  Foam  Cart  Data  Time  Range  Velocity,  No  Whitening  or  Extrapolation 
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Figure  46  Akela  Foam  Cart  Data  Time  Range  Velocity  After  Whitening  or  Extrapolation 
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During  a  collection,  individual  scatterers  on  the  boat  will  come  and  go  as  they  arc  shadow  ed/exposed,  and  their 
returns  wax  and  wane  due  to  angular  response.  The  algorithm  used  to  find  scattercr  paths  must  account  for 
these  effects.  We  chose  to  break  the  time  on  target  into  1 1  overlapping  regions.  Each  region  overlapped  half  of 
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the  previous  region.  Path  segments  were  found  within  each  region,  and  were  combined  and  ranked  based  on  the 
minimum  value  along  the  length  of  the  path  (higher  minimum  value  implying  a  better  segment).  Figure  47 
shows  the  scattcrcr  path  segments  that  were  extracted  on  top  of  a  projection  of  the  range  compressed  view  of 
the  time-range-velocity  cuboid. 


Figure  47  Akela  Foam  Cart  Data  Scatter er  Paths  Shown  on  Range  Compressed  signals 


Theoretical  Considerations 

Baseline  scenario 

Consider  a  sensor  at  coordinates  (0,0, h).  Consider  a  three-dimensional  target  whose  center  travels 
along  a  straight  path  defined  by  (x  +  vt,y,z),  wherein  t  denotes  time  and  v  denotes  velocity.  This 
situation  is  illustrated,  from  various  points  of  view,  in  the  follow  ing  figure,  for  y=10  and  z=10. 
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Various  views  of  the  sensor,  target,  target  path,  and  rays  from  sensor  to  the  target  path. 

The  rays  from  the  sensor  to  the  target  path  all  lie  in  a  fixed  plane,  the  plane  determined  by  the  sensor 
location  and  the  line  along  which  the  center  of  the  target  proceeds.  In  a  far-field  scenario  (whereby 
we  mean  that  wavefront  curvature  is  negligible  on  the  scale  of  the  target  dimensions),  this  plane 
constitutes  the  slant  plane,  and  only  one  specific  two-dimensional  image  can  be  produced  That  is  the 
standard  range,  cross-range  SAR  (or  1SAR,  as  the  distinction  depends  only  on  ones  frame  of 
reference)  image,  in  which  the  direct  reflections  from  all  scattering  centers  are  orthogonally  projected 
into  the  slant  plane  Under  the  far  field  approximation  (planar  wave-fronts,  on  the  scale  of  the  target, 
from  each  pulse),  all  points  on  the  target  that  project  onto  the  same  point  in  the  slant  plane,  have  the 
same  far  field  range  histories  throughout  the  aperture;  hence  they  contribute  indistinguishable  signals 
to  the  radar  data. 
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“Apparent  roll” 


In  the  near  field,  both  the  aspect  and  elevation  angles  of  the  rays  from  the  target  center  to  the  sensor 
(and,  isomorphically,  the  aspect  and  depression  angles  of  the  rays  from  the  sensor  to  the  target  center) 
change  significantly  as  the  target  moves  along  its  path.  For  the  given  scenario,  these  angles  are 
plotted  in  the  following  figure: 

Aspect  and  eleva:  on  to  sersor  as  seen  by  tbe  boat 


Aspect  and  elevation  angles  to  the  sensor  as  seen  from  a  target  fixed  coordinate  frame. 

The  change  in  elevation  angle  gives  rise  to  the  beguiling  notion  of  “apparent  roll”  of  the  boat  relative 
to  the  sensor.  This  “apparent  roll”  exists  even  though  there  is  no  roll  of  the  boat  relative  to  the  planar 
surface  of  the  wrater  Hence,  there  will  exist  in  the  radar  data,  some  information  about  the  height 
(vertical  distance  from  the  center  of  the  target)  of  any  scattering  center  on  the  target. 

Ambiguity  in  three  dimensions 

The  “apparent  roll”  observation  seems  to  suggest  that  there  is  more  three-dimensional  information  in 
the  near  field  radar  data,  concerning  a  linearly  moving  target,  than  there  is  in  the  far  field 
approximation  of  the  measurements,  wherein  only  two-dimensional  measurements  (in  the  slant  plane) 
can  be  made  (given  only  the  information  in  the  data).  But  this  conclusion  is  troublesome  from  a 
topological  perspective. 

The  transformation  from  the  far  field  case  to  the  near  field  case  involves  a  smooth  deformation  of  the 
metric  on  the  space,  as  long  as  the  target  stays  away  from  the  only  singularity,  which  is  at  the  sensor. 
So  as  long  as  the  target  does  not  collide  with  the  sensor,  topological  invariants,  such  as  the  dimension 
of  the  slant  plane,  and  that  of  its  complement  (the  span  of  its  normal  vector)  should  be  the  same  as  in 
the  far  field  case.  So  the  dimensionality  of  the  resolvable  target  shape  information  should  be 
conserved. 

So  where  does  the  ambiguity  (represented  in  the  far  field  by  the  projection  along  the  normal  to  the 
slant  plane)  go?  The  answer  is,  it  bends  around  a  line  through  the  sensor  location,  which  is  parallel  to 
the  target  trajectory.  To  illustrate  this,  we  can  examine  the  ranges  (exact  near  field  ranges,  no 
approximations)  to  some  distinct  points  along  the  circular  arc  shown  in  red  in  the  following  figure. 
Each  of  the  seven  marked  points  along  the  red  curve  is  1  meter  from  its  nearest  neighbor. 
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A  rigid  arc  of  ambiguous  possible  scattering  center  locations  attached  to  the  rigid  boat  model. 
The  ranges  to  seven  points  along  this  arc  are  plotted  in  the  following  figure. 


Range  histories  from  7  distinct  points  along  an  arc  of  ambiguity  within  the  target  model  space.  All  7 
curves  exactly  coincide  throughout  the  data  collection. 

There  is  only  one  curve  apparent  in  the  figure,  because  all  seven  curves  are  exactly  the  same  Hence, 
the  radar  data  will  not  contain  information  with  which  to  distinguish  the  locations  of  the  seven 
scattering  centers. 


More  generally,  the  complete  set  of  arcs,  through  the  target  space,  all  of  which  are  parallel  to  the 
illustrated  arc,  forms  a  foliation  of  the  target  space,  in  which  leaf  of  the  foliation  (any  individual  arc) 
contains  a  continuum  of  possible  scattering  center  locations,  any  of  which  will  contribute  exactly  the 
same  signal  to  the  radar  data.  Hence,  the  complete  set  of  available  information  about  the  size  and 
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shape  target  is  still  two-dimensional,  not  three-dimensional.  In  fact  it  can  still  be  represented  as  an 
image  in  a  two-dimensional  slant  plane.  The  only  difference  from  the  far  field  approximation  case  is 
that,  in  the  near  field,  three-dimensional  scattering  center  locations  “project”  into  the  slant  plane  by 
sliding  along  the  arcs  of  the  foliation  rather  than  along  the  normal  vectors  to  the  slant  plane.  Even  in 
the  near  field,  to  be  able  to  put  unambiguous  three-dimensional  target  size  and  shape  information  into 
monostatic  radar  data,  the  target  must  deviate  from  a  straight  line  path. 

The  orthographic  projections  in  the  following  image  illustrate  an  example  of  the  ambiguity.  The  set  of 
blue  scattering  center  locations  are  all  on  the  surfaces  of  the  target  model.  But  reflectors  at  the  set  of 
red  scattering  center  locations,  traveling  as  a  rigid  body  along  the  same  trajectory,  would  generate 
exactly  the  same  radar  signal  history.  The  red  set  was  chosen  by  choosing  a  point  at  random  from 
each  leaf  of  the  foliation  which  contains  a  blue  point. 


X  ’ 


.  y 


Orthographic  views  of  scattering  center  locations  attached  to  the  target  (blue)  and  an  alternative  set  of 
rigid  scattering  center  locations  (red),  which  would  reflect  the  same  radar  data  set  from  a  linearly 
moving  target. 

For  an  algebraic  proof,  recall  that  h  is  the  height  of  the  sensor.  Recall  that  each  scattering  center 
travels  from  some  starting  position,  (x,  y,  z),  along  the  path,  (x(t),  y(t),  z(t))  =  (x  +  vt,  y,  z).  For  any 

particular  scattering  center  starting  position,  (x1#y1#  zx),  let  R1  =  yjy \  +  (za  —  h)2).  Then  let 
(x2,  Y2>  zt)  —  (xi>  Ri»  cos(0) ,  h  +  Rv  sin(0)),  for  any  arbitrary  value,  0.  Under  the  prescribed  rigid 
target  motion,  the  first  point  will  traverse  the  linear  trajectory,  (x1  4-  vt,y1#  Zj),  and  the  second  point 
will  traverse  the  linear  trajectory,  (x2  +  vt,y2,  z2).  At  any  time,  t,  the  range  from  the  sensor,  at 
location,  (0,0, h),  to  the  first  scattering  center  location  will  be 

ri  (t)  =  Vxi(t)2  +  yi(t)2  +  (z,(t)  -  h)2 

=  VOi  +v02  +y?  +  (zi  -h)2 

and  the  range  from  the  sensor  to  the  second  scattering  center  location  will  be 
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r2(t)  =  V x2  (t)2  +  y2  (t)2  +  (z2  (t)  -  h)2 


(x2  +  vt)2  +  yl  +  (z2  -  h)2 


=  7 (xj  +  vt)2  +  (R,  cos(0))2  +  ((h  +  R,  sin(0))  —  h)2 
=  7 (xi  +  vt)2  +  (Rj  cos(0))2  +  ((Rj  sin(0)))2 


(Xj+vt)2  +  Rj(cos(0)2  +  sin  (0)2) 


(x,  +  vt)2  +  R2 


(xa+vt)2  +  (Jyf  +  (zj  -  h)2)2 


(xj  +  vt)2  +y\  +  (Zi  -  h)2 


=  V*i(t)2+yi(t)2  +  zj(t)  -  h)2 

=  ri(t) 

So  for  any  value,  0,  and  any  time,  t,  the  signal  from  the  second  scattering  center  position  will  be  the 
same  as  the  signal  from  the  first  scattering  center  position 


The  fact  that  the  linear  motion  scenario  generates  a  signal  data  set  that  does  not  support  three- 
dimensional  imaging  really  should  not  seem  surprising.  How  else,  given  data  generated  from  such  a 
situation,  could  we  have  used  the  two-dimensional  migration  program  to  form  a  perfectly  focused 
two-dimensional  image  of  the  simulated  data  (which  was  made  from  a  fully  three-dimensional  target 
model)?  If  there  were  genuine  three-dimensional  scattering  center  location  information  in  the 
simulated  data,  the  two-dimensional  migration  algorithm  would  be  unable  to  fully  focus  all  the  signal 
energy. 

The  target's  perspective 

It  can  be  easier  for  to  analyze  some  aspects  of  the  situation  from  the  inertial  reference  frame  of  the 
target.  Imagine  an  observer  inside  the  target,  who  does  not  see  out,  but  it  who  is  able  to  continuously 
measure  the  direction  and  distance  from  which  the  radiation  is  arriving.  This  observer  can  only 
quantify  these  directions  of  radiation  arrival  w  ith  respect  to  some  reference  frame  w'hich  is  defined  by 
objects  internal  to  the  target.  From  the  observers  world  view,  he  and  his  reference  frame  are 
stationary,  and  changes  in  the  direction  of  arrival  of  the  radiation  must  be  attributed  to  movement  of 
the  external  source  of  the  radiation. 


For  our  scenario  the  observer  will  find  that  the  aspect,  elevation,  and  distance  from  which  the 
radiation  comes  all  change  over  time.  But  despite  those  changes  in  aspect  and  elevation,  all  the 
illumination  directions  will  fit  in  a  common  plane,  containing  the  origin  defined  by  the  internal 
observer. 


6.  Description  of  the  McCorkle  Antenna 

The  antenna  that  was  designed  on  this  program  was  intended  to  demonstrate  that  designs  exists  that  could  fit  in 
very  small  form-factors,  be  conformal,  and  still  have  the  gain  and  other  characteristics  needed  for  this  type  of 
application.  A  printed  antenna  of  a  cardioid  type  was  custom-designed  by  John  McCorkle  for  this  program.  This 
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type  of  antenna  can  be  both  broadband,  and  have  reasonable  flat  gain  over  that  band.  The  physical  antenna 
produced  is  shown  in  Figure  48  below. 


Figure  48:  Antenna  pair  of  cardioid  type  produced  by  John 

McCorkle. 


This  shows  two  such  antennae  mounted  on  a  backplane.  The  overall  package  measures  1  r’x7”x3”.  The  antenna 
frequency  dependent  gain  and  voltage  standing  wave  ratio  (VSWR)  performance  is  shown  in  Figure  49  and 
Figure  50  below. 


w 


* 

i 


Figure  49:  Frequency  dependent  gain  plot  of  McCorkle  antenna  pair . 


It  can  be  seen  that  the  gain  of  this  antenna  in  the  band  above  800MHz  is  6dB.  The  VSWR  ranges  from  1 .25  to 
2.5  over  this  same  band.  This  means  that  the  efficiency  of  this  antenna  is  relatively  low,  but  is  generally 
acceptable.  Because  of  the  poor  match  at  some  frequencies,  it  is  important  to  minimize  cable  lengths  to  avoid 
having  reflected  power  appear  in  ranges  of  interest.  This  type  of  antenna  can  be  tuned  to  be  very  efficient  in 
bands 
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7.  Analysis  of  Akela  RADAR  for  Collection  Suitability 

The  purpose  of  this  section  is  to  document  the  procedure  used  to  establish  the  signal-to-noise  ratio  (SNR)  and 
noise  power  for  the  Akela2  radar  unit  The  lack  of  information  regarding  the  noise  figure,  or  knowledge  of  the 
internal  structure  of  the  receiver  means  a  straight  forward  derivation  from  the  radar  equation  is  not  available 
For  this  reason,  measurements  were  made  of  a  target  w  ith  known  radar  cross  section  (RCS),  at  a  known  range, 
and  the  system  noise  power  was  estimated  using  the  radar  range  equation.  Using  this  established  noise  power, 
SNR  curves  vs.  range  were  generated  for  targets  of  various  cross  section  in  order  to  ensure  that  adequate  SNR 
could  be  realized  in  upcoming  data  collection  experiments  given  the  required  ranges  and  RCS. 

Signal-to-Noise  Ratio 

To  make  reliable  radar  measurements,  the  received  signal  power  must  be  significantly  above  the 
receiver  noise.  This  is  usually  quantified  in  terms  of  a  signal-to-noise  ratio  The  standard  equation  used  for 
monostatie  radar  measurements  is: 


SNR= 


pT& 


m  ^  P  vj_  Ppa  „  E3-y»  □ 

4nR2TLT  47iR^  4nkBT  iysBnFLs 


where 


Pt 

average  transmit  power 

r 

transmit  pulse  width 

PRI 

pulse  repetition  interval 

Rt 

range  to  target 

Gt 

transmit  antenna  gain 

Gr 

receiver  antenna  gain 

CT 

radar  cross  section 
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Lt  transmitter  loss 

F  receiver  noise  figure 

Ls  system  losses 

A  wavelength 

Gsp  processing  gain 

kB  boltzman’s  constant 

Tsvs  system  temperature 

Bn  noise  bandwidth 

This  can  be  generalized  into  the  form: 

P 


SNR- 


where, 


□jiCKJ  P„ 


Ps  received  signal  power 

Pn  receiver  noise  power 

and, 

F  =  rTtipRWTGR^GslLs 

^n=^TnsB„F 

By  simply  substituting  the  appropriate  parameters  into  the  above  equation,  the  Signal-to-Noise  ratio 
(SNR)  can  be  predicted  to  an  accuracy  commensurate  with  the  combined  aeeuraey  of  the  parameters  entered. 
Unfortunately,  the  internal  details  of  the  Akela2  radar  hardware  are  unknown,  so  the  receiver  noise  figure  and 
noise  bandwidth  are  not  known. 

However,  by  making  a  simple,  monostatie  measurement  of  a  target  of  known  radar  cross  section,  at  a 

P 

known  range,  with  antennas  of  known  gain,  and  using  a  known  processing  gain,  the  system  noise  power  (  ”  ), 

k  T  B  FL 

consisting  of  the  B  s>'s  n  S  term  in  the  equation,  ean  be  calculated  from  the  measurement  results. 

To  begin  with,  it  is  important  to  make  a  worst  ease  estimate  in  order  to  establish  bounds  on  the 
operational  range  for  targets  of  various  expected  radar  cross  sections.  F  urther,  it  is  important,  in  the  near  term, 
to  be  able  to  estimate  this  performance  for  the  same  conditions  under  which  the  radar  will  operate,  meaning  that 
the  same  parameters,  such  as  transmit  and  receive  gate  timing,  should  be  used  in  the  measurement  to  establish 
the  noise  power,  and  hence  included  in  the  parameters  of  the  model.  To  do  this,  the  generalized  equation  above 
ean  be  simplified  somewhat. 

First,  the  angle  dependent  antenna  gain  terms  ean  be  replaced  by  constant  values  that  reflect  the  worst 
ease  expected  gain  for  both  the  transmit  and  the  receive  antennas.  In  the  noise  power  measurement  the  target 
will  be  located  on  the  boresight  of  the  antennas.  Further,  the  antenna  beams  to  be  used  in  the  data  collections 
have  extremely  broad  beams,  and  will  thus  have  little  variability  in  their  gain  patterns  within  the  scene  For 
upcoming  experiments,  the  Large  MeCorkle  antennas  will  be  employed,  and  have  a  peak  gain  of  6  dBi.  In  order 
to  consider  the  worst  ease  conditions  in  generating  the  predicted  SNR  performance,  the  -1  dB  point  of  the 
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pattern  (5  dBi)  will  be  used  as  the  intended  collection  scenarios  will  likely  limit  secnc  coverage  to  this  area. 

Both  antennas  will  be  the  same,  so  this  value  will  be  used  for  both  transmit  and  receive  antenna  gain. 

The  transmitter  loss  term,  Lr,  will  be  combined  into  the  average  transmit  power  term,  Pr,  with  the 
understanding  that  any  losses  due  to  cabling,  connectors,  etc.  between  the  radar  output  and  the  transmit  antenna 
must  be  reflected  in  the  Pr  term  in  the  equation.  The  Akela2  is  capable  of  putting  out  50mW  CW  power 
(nominally).  To  account  for  the  cables  used  in  the  experiments,  I  dB  of  loss  will  be  included.  Also,  in  order  to 
account  for  some  variability  over  the  band,  I  dB  of  loss  will  be  included  so  that  the  Prtcrm  will  be  31 .5  mW 
(15  dBm). 

Also,  since  the  measurements  are  intended  to  be  made  using  the  transmit  and  reecive  gating,  the  power 
will  be  reduced  by  the  duty  faetor.  For  these  experiments,  the  transmit  “on”  time  is  30  nsec,  while  the  time 
between  pulses,  or  pulse  repetition  interval  (PR1)  is  171  nsee. 

The  wavelength  will  be  estimated  from  the  middle  of  the  frequency  band  Since  for  the  expected  data 
collections  the  band  is  expected  to  cover  500  to  3000  MHz,  the  wavelength  used  here  will  be  17.1  em. 

The  processing  gain  in  the  equation  consists  of  the  usual  time  bandwidth  produet  realized  from  the 
range  compression  process,  w  hich  is  the  number  of  FFT  points  processed  (also  the  number  of  frequency  steps  in 
the  waveform).  The  desired  transmitted  bandwidth  as  stated  above  is  2500  MHz,  and  for  the  calibration 
measurements  in  the  chamber,  the  ambiguous  range  is  set  to  30  meters  In  the  field  experiments,  the  ambiguous 
range  will  be  set  to  100  meters.  This  means  that  frequency  step  size  in  the  stepped  linear  FM  radar  signal  for  the 
chamber  measurements  is  5  MHz,  resulting  in  501  steps/'FFT  points.  For  the  field  experiment,  the  frequency 
step  size  will  be  3  MHz  and  there  will  834  steps,  hence  834  FFT  points  for  the  range  compression. 

While  a  target  in  the  chamber  can  be  well  controlled  as  to  its  range  placement,  in  the  field  experiments 
the  targets  may  be  loeated  proximate  to  a  range  resolution  cell  boundary,  and  may  therefore  spread  the  target 
energy  over  2  range  resolution  cells.  To  compensate  for  these  effects,  a  faetor  of  2  processing  gain  reduction  is 
introduced  to  aceount  for  the  spreading  across  range  resolution  cells  in  the  SNR  prediction  for  the  field 
experiments,  or  41 7. 

System  receiver  losses  are  estimated  to  be  restricted  to  cable  losses,  and  at  these  low  frequencies,  are 
expected  to  be  limited  to  I  dB. 

To  estimate  the  noise  power  for  the  Akcla2  radar,  a  measurement  was  taken  of  a  sphere  with  RCS 
1 7.8  dBsm  at  a  range  of  6. 1  meters  in  the  MTRI  antenna  ehamber.  Clutter  in  the  chamber  environment  due  to 
the  imperfect  absorber  at  the  low  frequency  of  the  operation,  and  primarily  due  to  returns  from  the  metal  ceiling 
of  the  chamber  which  cannot  be  covered  w  ith  absorber  due  to  fire  code  restrictions,  as  well  as  the  internal  noise 
generated  from  feed-through  of  the  transmitted  signal  internally  and  at  the  antennas,  requires  a  background 
subtraction  proeess  in  order  to  enhance  the  target  return.  Another  measurement  of  the  ehamber  w  ith  sphere 
removed  was  taken,  and  coherently  subtracted  from  the  sphere  measurement.  The  resulting  range  response 
(averaged  over  all  the  sweeps  taken)  is  shown  in 
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Ave  Data  Rng  Profile  -  Ave  BG  Rng  Profile  from  Rotary  Platform  with_sphere_lower_forbg4_dat< 
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Figure 


Measuring  from  the  peak  of  the  target  at  6. 1  meters  to  the  nearby  surrounding  unoiscM  resulted  in  an 
estimated  SNR  for  this  target  of  14  dB.  The  measurement  was  done  by  eye,  and  eould  possibly  be  improved  by 
integrating  the  noise  over  several  range  bins  on  either  side  of  the  target  return  (this  wil  l  be  considered  for  future 
work).  This  number  was  plugged  into  the  rearranged  range  equation  to  estimate  the  system  noise  power  The 
parameters  used  for  this  estimate  are  summarized  below: 


PT  =  31.5  mW  average  transmit  power 

1  ~  30  nsee  transmit  pulse  width 

PR  I 

=  171  nsee  pulse  repetition  interval 

Rt  -  6. 1  m  range  to  target 

Gt  5  dBi  transmit  antenna  gain 

G/t  =  5  dBi  receiver  antenna  gain 

cr  =  -1 7.8  dBsm  radar  cross  section 
Ls  -  1  dB  system  losses 

A  =  17.1  em  wavelength 

Gy)  =  510  processing  gain 

And  the  resulting  “measured'1  noise  power  for  the  system  is  ealeulated  to  be: 
P"  =  3.0968  X  I  O'6  mW. 
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Figure  51:  Range  Profile  of -17.8  dBsm  Sphere  with  (blue  curve)  and  without  (red  curve) 

chamber  background  subtracted 


Using  the  noise  power  above  and  the  parameters  for  the  expected  field  collections  as  specified  below 
(discussed  above),  a  set  of  SNR  vs.  range  curves  for  several  different  RCS  targets  were  generated  and  arc 
shown  in 
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Akela  Radar  with  500  mW  Xmit  Power  Amp 
SNR  vs  Range  for  RCS=-20  to  20  dBsm 


Figure  52.  Note  that  transmit  power  reflects  the  use  of  a  power  amplifier  that  boosts  the  output  power  to  a  peak 
of  500mW  that  will  be  employed  in  the  field  measurements. 
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Akela  Radar  with  500  mWXmit  Power  Amp 
SNR  vs.  Range  for  RCS=-20  to  20  dBsm 
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Figure  52.  Predicted  SNR  vs.  Range  Performance  for  the  Akela  Radar  with  500  mW  Power 
Amp  Based  on  Noise  Power  Calculated  from  Sphere  to  Chamber  Noise  Level 


Note  that  the  “measured'1  SNR  that  resulted  in  the  noise  power  measure  above  is  actually  based  on  a 
signal-to-clutter  level,  and  may  not  accurately  depict  expected  performance  under  field  conditions  where  the 
clutter  is  far  more  benign  than  that  of  the  chamber.  A  different,  potentially  more  accurate,  measure  of  SNR  was 
sought.  The  radar  was  operated  with  the  same  parameters  as  for  the  calibrated  sphere  measurement  above,  with 
the  transmit  and  receive  ports  terminated  by  50Q  loads.  The  resulting  data  were  averaged  over  each  frequency 
bin,  and  this  average  was  coherently  subtracted  from  the  data  to  remove  the  feed-through  signal,  yielding  a 
better  estimate  of  the  noise  level.  Since  the  radar  parameters  were  the  same,  the  resulting  range-compressed 
outputs  arc  on  the  same  scale,  and  the  SNR  of  the  sphere  can  be  estimated  by  measuring  the  noise  level  from 
this  curve  ( 
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Residual  Range  Profile  from  Antenna  Chamber  int_noise_port_1T_2R_10_28_2010_data 
Range  Profile  -  Average  Range  Profile 


Figure  53), 
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Residual  Range  Profile  from  Antenna  Chamber  int_noise_port_1T_2R_10_28_2010_data 
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Figure  53:  Compressed  Range  Profiles  with  Average  Range  Profile  Coherently  Subtracted 
for  Aklea  Radar  with  All  Ports  Terminated 
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Figure  53,  the  SNR  of  the  calibrated  sphere  target  is  estimated  to  be  2ft  dB.  Using  this  to  calculate  the  noise 

power  yields  ^  n  —  1 .2329  X  10  mW,  and  with  the  same  parameters  as  above,  the  SNR  vs.  range  curves  arc 
re-calculated  and  shown  in  Figure  54  below. 
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Akela  Radar  with  500  mW  Power  Amp 
SNR  vs  Range  for  RCS=-20  to  20  dBsm 


Figure  54:  Predicted  SNR  vs,  Range  Performance  for  the  Akela  Radar  with  500  mW  Power 
Amp  Based  on  Noise  Power  Calculated  from  Sphere  to  Terminated  Radar  Noise  Level 


Recently,  it  was  found  that  there  is  a  non-linear  relationship  between  the  Akela  radar’s  noise  figure  and 
the  receiver  gain  value.  This  gain  is  in  an  amplifier  downstream  from  the  LNA.  The  following  statement  is  from 
Allan  Hunt  from  Akela;  uFor  the  ease  where  you  are  using  the  pulse  modulator  with  the  LNA,  the  SNR 
improves  as  you  move  from  a  gain  setting  of  0  to  6  with  the  improvement  being  about  10  dB  at  gain  6.  Above 
that,  the  SNR  starts  to  decrease.” 

Unfortunately,  the  above  model  and  measurements  assumed  a  more  well  behaved  linear  relationship  in 
the  receiver  with  any  gain  being  equally  applied  to  the  signal  and  noise.  This  is  clearly  erroneous,  and  based  on 
the  above  statement,  and  that  the  measurements  here  were  taken  with  gain  values  of  8  (=  24  dB),  it  is  estimated 
that  the  noise  power  calculated  above  will  be  reduced  by  5  to  6  dB,  and  the  SNR  curves  would  be  improved  by 
the  same  amount.  It  is  clear  that  more  measurements  will  be  required  to  get  a  more  accurate  model  for  the  Akela 
SNR  These  will  be  undertaken  as  eost/sehedule  restrictions  permit 
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8.  Summary 

This  report  outlines  the  progress  made  through  December  2010  of  the  program.  During  that  reporting 
period,  simulations  and  local  collections  occurred.  Both  the  simulations  and  local  data  collections 
provided  useful  data  for  the  algonthm  development  tasks.  In  addition,  a  new  antenna  system  was 
designed,  tested,  and  fabricated  for  use  in  the  field  collections  which  will  occur  in  January,  2011. 

Progress  was  made  on  the  signal  processing  algorithm.  3D  reconstructions  of  simulated  and  locally 
collected  data  were  created,  In  the  instance  of  far  field  scenarios,  the  3D  scattering  center  estimates 
were  quite  good. 

Our  plan  is  to  collect  true  river  boat  data  in  the  coming  period  and  perform  3D  reconstruction  on  that 
data.  The  data  collection  will  occur  in  January  201 1  on  the  Pearl  river  in  southern  Mississippi 
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Vol.  II:  Information-Theoretic  Analysis  &  Performance  Bounds 
for  Super-  Resolution  (SR)  Video  Imagery  Reconstruction 


1  Introduction 


Video  and  wide-area  surveillance  imaging  systems  acquire  massive  amounts  of  data  owing  to 
demands  for  high-frame  rates  and  large  fields-of-view.  Collected  data  must  he  transferred  to 
analysts  through  channels  with  finite  transfer  rate  capacities.  The  large  mismatch  between  the 
volume  of  the  data  collected  and  the  rate  at  which  it  can  be  transferred  results  in  long  delays  between 
data  acquisition  and  data  analysis.  We  are  proposing  a  joint  collection-processing  method  where  the 
volume  of  collected  data  can  be  reduced  while  maintaining  large  fields-of-view  (FOV)  and  high 
video  frame  rates  by  acquiring  fewer  pixels,  i.e.,  subsampling.  Subsampling  by  a  factor  of  2  or  4 
(equivalent  to  reducing  resolution  by  a  factor  of  2  or  4),  and  maintaining  the  same  tield-of-view, 
results  in  an  overall  reduction  in  data  volume  by  factors  of  4  and  16,  respectively.  Super-resolution 
image  reconstruction  algorithms  provide  a  means  of  recovering  fine  detail  from  low  resolution 
imagery,  particularly  video  imagery.  That  is,  in  principle,  data  volumes  can  be  reduced  without 
sacrificing  resolution.  Alternatively,  SR  image  reconstruction  also  provides  a  means  for  increasing 
imaging  system  FOVs  while  maintaining  resolution  and  current  data  sizes.  The  limits  of  line  detail 
recovery  from  lower-resolution  (aliased)  imagery  depend  on  various  factors: 

•  Scene  information  content 

•  Sub-sampling  factor  (relative  to  Nyquist) 

•  Optical  and  motion-induced  blur  (known  and  unknown) 

•  Signal-to-noise  ratio 

•  Camera  frame  rate 

In  the  project,  we  are  developing  an  information-theoretic  framework  for  super-resolution  that 
will  provide  quantitative  limits  of  achievable  super-resolution  image  reconstruction  under  general 
imaging  conditions  (including  ’’warping”  effects  of  a  changing  camera  perspective)  as  a  function  of 
key  variables. 

Our  approach  for  development  of  this  formal  theory  centers  on  deriving  appropriate  models  for  (i) 
estimation-theoretic  lower  bounds  (e.g.,  Cramer-Rao  bounds)  and  (ii)  Shannon  mutual-information 
between  fine-resolution  scenes  and  collected  data,  with  the  preliminary  emphasis  in  the  early  part 
of  the  project  being  the  former.  Both  of  these  provide  a  rigorous  framework  for  quantitatively 
assessing  performance  of  super-resolution  (SR)  video  imagery  reconstruction  methods  in  recovering 
the  desired  fine-resolution  detail.  Using  these  models,  we  will  be  deriving  bounds  (for  either 
estimation  accuracy  or  mutual  information)  as  a  function  of  key  scene  and  acquisition  parameters. 
The  bounds  will  can  incorporate  spatial  frequency  and  temporal  domain  weightings  that  emphasize 
the  particular  goal,  e.g.,  visual  fidelity,  or  registration,  or  subtle  change  detection.  The  main  focus 
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of  the  present  project  is  for  video  monochromatic  and/or  color  data  but  the  framework  is  flexible 
enough  to  encompass  hyper-spectral  and  polarimetric  imagery  as  well. 

The  organization  of  this  is  as  follows:  the  next  Section  2  provides  a  description  of  the  progress 
for  the  just-completed  quarter  of  the  project.  Section  3  provides  a  more  detailed  description  of  the 
mathematical  framework  and  the  associated  information-theoretic  Cranier-Rao  bounds.  Finally 
Section  4  provides  a  summary  of  accomplishments  and  plans  for  work  in  the  upcoming  quarter. 


2  Progress  in  the  Fourth  Quarter  of  2010 


A  high-level  list  accomplishments  for  the  quarter  just  completed  is  given  in  Table  1 . 


Information  content  metrics  for  video  imagery: 

1.  Developed  general  theory  for  Bayesian  Cramer-Rao  bounds 

2.  Incorporation  of  optical/detector  blur,  aliasing/sampling,  and  noise  models 

3.  Matlab  implementation  complete 

4.  Translation  only  case  -  “arbitrary”  image  sizes 

5.  General  case  -  smaller  images 

In  simulation  compared  theoretical  metrics  with  restoration  accuracy: 

1.  Translation-only  case/larger  images 

2.  Full  perspective  ease/smaller  images 

3.  Various  levels  of  sampling  and  SNR 

4.  Results  look  quite  promising 


Table  1:  Accomplishments  of  Video  Super-Resolution  Project  -  Quarter 

We  have  carried  out  mathematical  development  for  the  Bayesian  Cramer-Rao  bounds  (BCRB)  for 
video  imagery  -  the  full  mathematics  of  much  of  this  framework  is  included  in  Section  3  and  the 
appendices.  In  particular,  we  have  developed  a  mathematical  framework  for  deriving  estimation- 
theoretic  bounds  on  image  restoration  performance  as  a  function  of  number  of  frames,  aliasing, 
signal  levels,  and  camera  movements  during  the  collection,  and  thereby  provide  our  main  metric  for 
image  restoration  quality/information  The  mathematical  analysis  for  the  single  channel  is  complete 
and  documented  within  a  memo/paper  form,  and  a  significant  portion  of  this  is  included  in  the 
Section  3.  A  significant  feature  of  our  model  is  that  it  is  totally  general  with  a  forward  imaging 
model  that  includes  optical  blur,  detector  blur,  aliasing,  Poisson  shot  noise,  and  detector  noise.  We 
have  implemented  two  forms  of  the  single-channel  BCRBs  in  MATLAB.  One  is  an  image-domain 
method  and  allows  for  a  general  camera  perspective  -  the  limitation  of  this  approach  is  that  it  is 
only  applicable  to  smaller  100x100  chip  sizes  which  may  be  useful  for  smaller  regions  of  interest 
applicability.  We  have  also  implemented  a  frequency-domain  method  -  it  is  able  to  operate  on 
full  scenes,  with  the  scene  size  only  limitation  being  that  of  memory,  but  it  is  currently  limited  to 
only  translation  methods.  Current  research  is  focused  on  deriving  extension  of  these  frameworks 
to  the  more  general  case  of  large  scene  sizes  and  realistic  evolution  of  camera  perspectives  that 
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are  more  general  than  translations.  The  main  challenges  here  are  computational  and  that  is  the 
reason  for  the  significant  research  effort  in  this  regards.  We  have  derived  a  full  Bayesian  maximum 
a-posteriori  (MAP)  image  dealiasing  and  restoration  algorithms  for  multi-frame  video  -  the  output 
results  from  this  appear  to  be  quite  good.  The  current  version  of  this  algorithm  relies  on  knowledge 
of  the  camera  perspectives,  and  has  the  nice  feature  that  the  estimators  have  a  closed  form  and  do 
not  require  an  iterative  optimization.  It  is  important  to  note  that  these  models  could  be  appropriate 
in  very  well-controlled  scenarios  where  the  approximate  camera  perspectives  are  well-known.  As 
mentioned  above,  our  mathematical  models/framework  is  general  in  that  it  allows  for  known  or 
unknown  camera  perspectives,  but  we  have  not  implemented  the  general  unknown-perspective  case 
in  code  and  tested  it  out.  Figure  I  shows  an  an  example  of  a  de-aliased  restoration  -  on  the  left  are 
the  original  images  and  on  the  right  is  the  restoration. 


Figure  1 :  Example  restoration  using  the  frequency  domain  method.  Left:  Four  video  frames  sampled 
at  0.25  the  Nyquist  rate.  Right:  Restoration  using  the  derived  closed  form  estimator. 


Another  accomplishment  on  this  code  is  that  we  extended  the  regularization  components  of  the 
dealiasing  image  restoration  code  to  allow  for  a  general  regularization  matrix  (Bayesian  Prior)  on 
the  image  -  this  allows  us  to  regularize  to  the  content  of  the  imagery  and/or  to  the  goals  of  the 
restoration,  i.e.,  focused  more  on  high-frequency  information  (e.g..  edges)  as  opposed  to  lower 
frequency  information  (e.g.,  texture).  We  also  derived  a  computationally  efficient  method  for  finding 
the  appropriate  magnitude  of  the  regularization  coefficient  that  corresponds  to  generating  minimum 
mean-squared  error  restorations.  This  is  useful,  as  it  is  in  some  sense  the  “optimal”  choice  with 
respect  to  squared  error  despite  being  a  function  of  the  unknown.  We  have  also  begun  exploration  of 
a  blind  regularization  selection  approach,  and  demonstrated  its  efficacy  as  compared  to  other  current 
techniques  (e.g.,  generalized  cross-validation  -  GCV)  on  some  simpler  problems  and  are  working 
on  extending  to  our  more  general  situation.  This  development  is  very  important  since  the  current 
data-adaptive  techniques  for  selecting  the  regularization  in  single/multi-frame  image  restoration 
are  not  a  reliable,  and  this  is  especially  true  in  the  more  challenging  scenarios  of  doing  signilicant 
de-aliasing.  To  have  a  reliable  technique  in  this  area  would  be  a  big  enabler  for  doing  accurate  and 
informative  image  restoration. 


11-3 


The  Bayesian  Cramer-Rao  bound  analysis  is  designed  around  having  a  known  point-spread-function 
(PSF)  model  for  the  imaging  sensor.  This  is  a  reasonable  working  assumption  only  when  relatively 
accurate  estimates  of  the  PSF  can  be  obtained.  Along  these  lines,  we  generated  an  algorithm/ap¬ 
proach  for  doing  estimation  of  the  system  blur  utilizing  a  known  target  that  is  relatively  optimized 
for  this  application  In  many  systems  one  is  able  to  collect  against  a  known  target  before  deploying 
a  sensor  and  this  algorithm  aids  in  the  design  of  that  target  to  ensure  accurate  PSF  estimates. 
Our  procedure/framework  utilizes  a  class  of  targets  (2D  extensions  of  the  famous  Cantor  sets 
in  Mathematics)  and  these  were  chosen  so  as  to  provide  significant  robustness  to  mismatches  in 
the  assumed  sensor  model  (uncertainly  in  focal-length,  detector  pitch,  etc).  We  have  generated  a 
MATLAB  code  that  automatically  accounts  for  target  resolution  (and  includes  printer  capabilities) 
and  suggests  a  target  (both  its  content,  and  its  scale)  to  ensure  successful  experiments.  Once  you 
have  collected  data  off  of  this  target,  our  estimation  model  jointly  estimates  the  unknown  PSF 
subject  along  with  a  perspective  mismatch  that  will  always  be  there  and  needs  to  be  accounted  for. 
This  latter  feature  allows  for  a  user  to  not  require  that  the  target  be  perfectly  aligned  to  the  sensor's 
line-of-sight,  and  allows  for  reasonable  errors  in  the  assumed  sensor  parameters.  The  approach  was 
tested  and  verified  on  a  Sony  Power-Shot.  Figure  2  shows  a  version  of  an  original,  somew  hat  classi¬ 
cal  target  typically  utilized  (on  the  left)  and  on  the  right  is  the  new  target  which  we  propose.  Via 
simulations,  we  have  shown  our  target  is  a  significant  improvement  over  the  classical  for  PSF  estima¬ 
tion  and  it  is  also  much  more  robust  in  operations,  due  to  the  multi-resolution  properties  of  the  target. 


Figure  2:  (Left)  Classical  resolution  target  -  These  targets  tend  to  be  sensitive  to  collection  range 
due  to  their  concentration  of  spatial-frequency  energy  .  (Right)  Suggested  target  -  This  target  is 
designed  to  avoid  aliasing  in  the  printing  process,  but  has  asymptotically  uniform  spatial-frequency 
content 

Our  generalization  of  the  image  reconstruction  algorithms  and  information  metrics  to  imagery  taken 
at  an  arbitrary  sampling  rate  requires  an  “accurate”  generalization  of  the  aliasing  operator  (sampling 
of  collected  image  data  at  a  rate  that  is  sub-Nyquist)  and  the  “adjoint"  of  that  aliasing  operator, 
where  this  latter  operator  is  similar  to  the  adjoint/transpose  in  matrix  theory.  An  accurate  model 
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of  both  of  these  are  needed,  and  as  the  work  progressed,  we  recognized  that  our  current  models 
and  implementations  were  not  sufficiently  accurate  enough  for  the  demand  of  future  applications, 
especially  as  we  move  to  the  full-video  dealiasing  implementation.  So,  we  focused  on  developing, 
implementing,  and  verifying  aliasing  and  aliasing-adjoint  operators  that  can  generate  images  ac¬ 
quired  at  any  sampling  rate  Q  from  an  image  acquired  at  sampling  rate  P,  with  the  only  constraint 
being  2  >  P  >  Q.  As  mentioned  earlier,  the  details  in  developing  this  operator  are  quite  subtle,  and 
previous  implementations  had  errors,  some  significant  and  some  not  so  significant.  The  physical 
mechanism  that  is  being  simulated  can  be  described  in  the  following  way.  Given  an  image  acquired 
with  a  detector  array  consisting  of  N  x  M  pixels,  each  of  size  R  x  S  microns,  with  sampling  rate  PI ; 
generate  the  image  that  would  have  been  acquired  by  a  detector  of  sampling  rate  P2  constrained  to 
the  same  physical  array  dimensions  of  NR  x  MS  microns.  We  accomplished  this  by  analyzing  the 
aliasing  process  in  the  Fourier  domain.  In  the  Fourier  domain,  the  equivalent  spectrum  of  the  image 
sampled  at  rate  P2  is  obtained  by  appropriately  shifting  and  summing  the  frequency  components 
of  the  original  image  so  as  to  maintain  the  maximally  smooth  interpolated  signal  of  the  object. 
The  operator  was  generalized  to  include  (i)  resampling  with  and  without  aliasing  and  (ii)  energy 
preservation  under  both  the  LI  and  L2  norms.  In  addition  to  this  “forward”  aliasing  operator,  wc 
also  derived  and  implemented  the  full  adjoint.  The  adjoint  operator  is  quite  complicated,  but  it  is 
required  for  computation  of  the  information  metrics  as  well  as  gradient-based  optimization  routines 
for  image  reconstruction.  The  subtle,  but  important  changes  in  the  operators  are  essential  estimating 
“into  the  tail”  of  an  objects  spectral  content,  which  we  need  to  do  good  adaptive  generation  of 
information  quality  metrics  from  actual  collected  data.  In  Figure  3  is  an  example  of  the  data  as 
generated  with  this  new  code  at  5  different  sampling  factors,  with  Q=2  being  Nyquist  and  the  rest 
corresponding  to  different  undersanipling  factors. 


Figure  3:  Equivalent  images  sampled  at  (from  left  to  right)  Q=2,  1 .2,  0.85,  0.47,  and  0.27 

We  also  developed  a  mathematical  framework/approach  for  the  adaptation  of  the  current  approach 
to  information  content  metrics  for  color/multispectral  image  data.  The  main  focus  here  is  on 
the  transformed  space  of  luminance  and  chrominance.  We  do  not  believe  that  this  transforma¬ 
tion/framework  is  the  final  answer  in  the  complex  and  dynamic  world  of  color  imagery,  since 
there  are  competing  demands  for  “visual”  fidelity  as  seen  by  a  user  and  “quantitative”  fidelity  as 
demanded  by  an  exploitation  algorithm.  But  it  is  a  very  good  start,  and  it  is  expected  that  the 
“prior”  information  as  would  be  used  in  the  information  metrics  and  in  the  restoration  algorithms 
will  be  tailored  to  these  transformed  channels.  Again,  the  flexibility  of  our  framework  (arbitrary 
weightings  in  either  the  spectral  or  spatial  domain)  allows  us  the  ability  to  emphasize  cither  visual 
fidelity  or  exploitation-based  fidelity  within  the  video  color  image  restoration.  This  is  important 
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since  there  are  such  a  wide  and  varied  set  of  criteria  in  this  regards.  We  started  the  mathematical 
modeling  of  the  information  for  the  full  with  two  preliminary  concepts  having  been  identified  - 
the  first  is  based  on  the  use  of  sparse  matrix  representation  of  the  warp  operator  and  the  second 
is  based  on  doing  a  stitching  together  (essentially  aggregation)  of  subimage  information  metrics. 
We  are  in  the  process  of  fully  writing  out  the  mathematical  framework  of  both,  and  starting  the 
implementation  should  begin  next  month.  Also,  we  improved  the  accuracy  of  our  models  and 
simulations  by  the  incorporation  of  broadband  PSF  into  forward  and  inverse  model.  Real  gray-scale 
imagery  is  composed  of  light  from  across  the  visible  spectrum,  so  we  added  more  accuracy  (i.e., 
realism)  to  our  forward  and  inverse  models  by  developing  tools  that  can  automatically  generate 
broadband  point-spread  functions  consisting  of  an  arbitrary  weighting  of  spectral  components.  This 
tool  can  now  be  adapted  to  match  any  single-channel  image  (including  component  channels  of  a 
multispectral  image).  Figure  4  below  shows  a  version  of  two  point-spread-functions  (PSFs),  with 
the  one  on  the  left  corresponding  to  a  broadband  PSF  over  500-700  nm  range  and  the  one  on  the 
right  corresponding  to  a  narrowband  PSF  of  600  nm. 
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Figure  4:  Left  -  PSF  created  using  100  spectral  components  spanning  500-700  nm.  Right  - 
Monochromatic  PSF  at  600nm.  Both  PSFs  were  created  with  0.2  waves  of  RMS  wavefront  error  at 
600nm. 

We  have  carried  out  two  different  models,  one  in  the  so-called  Fourier  domain  and  one  in  the  image 
(spatial)  domain,  and  did  a  formal  and  rigorous  verification  that  these  two  methods  are  equivalent 
for  the  case  rigid  translations  for  a  wide  variety  of  targets  and  scenarios.  The  verification  of  the 
approximate  equivalence  of  the  two  methods  was  an  important  consistency  check  that  needed  to  be 
done  and  it  was  found  that  the  information  metrics  and  image  reconstruction  results  were  ~  3%, 
relative  difference  .  We  did  not  expect  the  two  methods  to  yield  identical  answers  due  to  numerical 
round-off  differences  between  the  methods  and  this  level  of  error  is  quite  acceptable  for  this  stage 
of  the  research,  and  moving  forward  with  our  algorithm  development  efforts. 
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3  Mathematical  Framework 


In  this  section,  we  provide  some  of  the  background  mathematics  that  are  critical  to  modeling  the 
forward  propagation  model  of  video  data,  and  the  inverse  estimation  framework  that  we  have 
developed  based  on  this. 


3.1  Forward  Model 

To  begin  the  section,  we  first  define  a  number  of  symbols  utilized  throughout  this  section. 


Definition  of  symbols 

f  =  Nyquist-samplcd  object  (Nf  x  Mf) 

W  =  geometric  warp  operator 
13  =  space-invariant  blur  operator 
A  =  aliasing  operator 


dj  =  jth  video  frame  (Nj  x  Md) 

Vj  =  geometric  warp  parameters  for  jth  frame 
G  f  =  guardband1  on  a  grid  of  size  (A/  x  Mj) 
G„  =  guardband2  on  a  grid  of  size  (Ar,/  x  M,,) 


Ej  =  noise  vector  for  the  jth  video  frame 

For  simplicity,  we  adopt  an  iid  Gaussian  noise  model,  the  jth  frame  in  the  video  sequence  is  written 


d  j  =  G(,AGfBW(vJ)f  +  £j.  (3.1 ) 

For  a  sequence  of  J  video  frames  we  have 

d  =  diag  ({GflAG/BW(vJ)}^1)  (l,  <g>  I(l/)  f  +  e.  (3.2) 

where  d  =  (dj , _ d j),  diag  {•} ' =1  is  a  block  diagonal  matrix  with  J  blocks,  -  (tq, . . . ,  vq), 

lj  =  ones(./.l),  tif  =  NfMf,  nd  =  NdMd,  I?1/  =  the  (n.f  x  n f)  identity  matrix,  and  e  = 
(ei . £j )• 


3.2  General  Bayesian  Cramer-Rao  Bound 


To  simplify  notation  we  define  Hv  =  diag  f  {GaAG/BW(vj)}j=1^  (lj  ®  Iri/)  where  v  =  (v| 
so  our  model  becomes 


d  =  Hvf  +  e. 


(3.3) 


For  simplicity,  we  assume  f  ~  A'(/tf,  £f)  where  Sf  is  positive  definite,  and  v  ~  N(fiv.  Sv)  where 
Ev  is  positive  definite.  Under  this  assumption  a  negative  log-likelihood  function  is 


where  '  indicates  conjugate  transpose  (or  Hermitian  transpose),  and  ||  •  ||  is  the  Euclidean  norm  on 

(C"rf. 


The  information  matrix  corresponding  to  (3.4)  is  a  sum  of  information  matrices  corresponding  to 
the  data-fit  term  and  the  prior  information  term,  S  r  =  So  +  Sp,  where  S y  is  the  total  information 
matrix,  Sp  is  the  information  matrix  corresponding  to  the  data-fit  term,  and  Sp  is  the  information 
matrix  corresponding  to  the  prior  term. 


The  information  matrix  corresponding  to  the  prior  term  is  straightforward  to  compute: 


V2  (^f'£f ’f+^v'E^v) 

®{nfxJq) 

\2  f  2  v  )\ 

P(Jqxnj)  Sv 

The  information  matrix  corresponding  to  the  data-ht  term  is: 


So  =  E 


d  -  Hv 


2 

A 

Sff 

Sfv' 

Svf 

Sw 

In  both  (3.5)  and  (3.6)  the  V2  operator  has  the  block  structure  shown  in  tigure  5. 


Up  Jq 


A 


np 

V2  = 

y 

A 

Jq 

V 


Vff 

<1 

Vjf 

V2 

v  vv 

(3.5) 


(3.6) 


Figure  5:  Block  structure  of  V2. 


We  are  ultimately  interested  in  the  upper-left  block  of  the  inverse  of  the  total  information  matrix, 
St  =  So  +  Sp.  By  applying  a  form  of  the  matrix  inversion  lemma  we  have 

[(Sff  +  Sf-1)-Sfv(Svv  +  Ev-1)-1Svf]-'.  (3.7) 


Next  we  apply  the  Woodbury  formula  to  obtain: 

=  [Sff  +  Sp1]-' 

object  independent 

+  [Sff  +  V] Sfv  ( [Sw  +  £-']  -  svf  [Sff  +  £-']  1  Sfv)  _1  Svf  [Sff  +  £(.  ']  ' 

^  v  ✓ 

object  dependent 


(3.8) 
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We  next  give  expression  for  Spr,  Sfv.  Svf,  and  Svv- 


S  a  =  E 


(3.9) 


Next  we  compute  the  columns  of  SfV.  Each  column  is  formed  by  first  computing  the  column 
gradient  of  L  with  respect  to  f,  and  then  performing  partial  differentiation  with  respect  to  the  warp 
parameters.  The  columns  are  arranged  according  to  figure  6.  So  the  vk  column  in  the  above  is 


4h'vhv 

a1 

jr2  (lS®In/)diag({W'(vj)B'G//A'G'(IGaAG/BW(vj)}^ ,)  (l,®I,t/) 


1 

Jq 

I 

1 

1 

T 

Sfv  = 

8v„(Vft)  8v„(V, /.)  . 

..  aVJj(vfL) 

_  *  * 

A 

Figure  6:  Format  of  Sfv  -  note  that  we  are  assuming  nj  »  Jq. 


given  by 
[Sfv]v^  =  E 


Avr^iid-Hvfir 


this  is  the  column  vector  corresponding  to  the  kth  warp  parameter  of  the  jth  frame  (3.10) 


=  E 


d 

dv 


jk 


[d  -  Hvf] 


^  ( -'-H'vE[d]  -  -H;Hvf  -  H'  tt— —  Hvf 
a2  \dvik  dvh  3vjk 


Tp~Hv  )  f. 
#vj* 


(3.11) 
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Now, 


diag  ({G„AGfBW(Vj)};=1)  (l,  <g>  In/) 


(3.12) 


=  (Hag 


(3.13) 


ro 


0 


g„aG/Bt^  w(vj)  ; 


(3.14) 


LO 


0 


So  we  have 


(3.15) 


Lastly,  noting  that  Svv  is  block  diagonal,  we  give  the  expression  for  the  jth  block  of  Svv' 


(3.16) 


(3.17) 


where  (•  •)  is  the  Euclidean  inner  product  on  C"d.  In  principle,  the  inverse  of  the  total  information 

matrix  can  now  be  computed.  While  the  diagonal  elements  of  the  inverse  total  information  matrix 
provide  mean-square  error  bounds  for  individual  parameters,  a  natural  global  scalar  performance 
measure  is  the  mean  of  the  diagonal  elements. 

3.2.1  Translations  only 

In  the  case  where  there  is  only  translation,  we  can  drop  the  guardband  operators  and  so  we  have  that 


In  the  case  where  there  is  only  translational  motion  between  the  video  frames,  our  performance 
bounds  can  be  computed  in  the  Fourier  domain.  This  can  be  done  via  the  following.  Let  Ta  denote 
the  unitary  version  of  the  2-D  FFT  operator  on  the  aliased  grid  and  let  T /  denote  the  unitary  version 
of  the  2-D  FFT  operator  on  the  fine  grid.  Note  that  T'a  and  T'j  are  the  normalized  inverse  2-D  FFT 
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operators.  For  this  case  we  use  the  following  notation  for  the  Fourier  counterparts  of  Sr-.  W.  B, 
and  A: 


Sr  =  j}SrJ> 

(3.19) 

W  = 

(3.20) 

B  =  Tf 

(3.21) 

A  =  TaATs 

(3.22) 

iNoting  that  since  B  is  block-circulant  with  circulant  blocks,  the  resulting  matrix  B  is  diagonal  and 
the  definitions  of  our  aliasing  operator  and  translation  operator  are  essentially  given  in  the  Fourier 
domain,  so  we  can  automatically  specify  the  exact  form  of  W  and  A.  Putting  this  with  the  obvious 
relation  of 


Fi?l  =  r,?r  =  h 

f.,K  =  K?*  =  i. 

where  we  note  that  the  '  denotes  the  conjugate  transpose.  Now  for  ,/  =  1,  easy  to  see  that 
Sfr  =  FjSffTf 

=  '(vi  )J>  FfB'TfFf  A'Fa  ra  ATfFfBTf  7}W(v,  )TS 

=  ^(^/W(v,)7»'  (FfB J»'  (rnATf)'(FaATf)  (T) BJ»  (^W(v,)J» 

=  -4w'(v,)B'A'ABW(v,)  (3.23) 

Note  the  warp  parameters  in  this  case  are  translation  parameters  and  there  are  two  per  video  frame. 
It  is  easy  to  generalize  the  above  to  the  general  case  of  J  >  2,  and  so  the  result  in  this  case  are  that 
the  elements  of  the  data-fit  information  matrix  become  (omitting  guardband  operations): 

Sr  =  ^  (1',  ®  I„,)  ding  ({wV;)B'A'ABW(Vj)}"=  J  (lj®In/),  (3.24) 

where  ~  denotes  the  Fourier  representation  of  the  quantities  and  operators\  We  note  the  operators 
<  W.  B  >  are  diagonal  and  that  the  product  A' A  can  be  permuted  into  block-diagonal  form. 


Also, 


[§f»k  -  4w'(Vj)B'A'AB^  W(v,)f 

=  -^W'(vJ)B'A'AB  (—2m  rampj  W(Vj)f 
[Sfv]vJy  =  ^yW'(vj)B'A'AB  (-2m  ramp„)  W(vj)f 

[Sw](vlx  VJy)  =  (AB(-2^'  rampj.)  W(vj)f  ,  AB  (~2m  rampy) 

'Fourier  representations  of  operators  are  obtained  by  pre  and  post  multiplication  by  a 
an  inverse  2D  DFT  matrix  respectively. 


(3.25) 

(3.26) 

(3.27) 

W(Vj)f) 
forward  2D 


<V  (3-28) 

DFT  matrix  and 


Note  that  [Svv]  is  block  diagonal  with  (2  x  2)  blocks. 


3.2.2  Performance  bounds  for  the  translation  only  case 


Recall,  the  total  inverse  information  matrix  consists  of  two  parts:  an  object  independent  (or  system) 
part,  and  a  part  that  depends  explicitly  on  the  object: 

[Sxlff1-  [Sfr  +  Sf'r1 

v  v  ^ 

object  independent 


+  [Sff  +  S^]  1  Sfv  ([Sw  +  £■']  -  Svf  [Sff  +  sp1]  1  Sfv)  'Svf  [Sff4-  Ef  '] 

V - v - 

object  dependent 

We  note  that  the  scale  of  the  object  is  incorporated  into  the  object-independent  part  through  the 
prior.  All  of  the  above  operators  have  Fourier  domain  representations,  and  there  is  an  overall  scale 

factor  of  —  from  the  discrete  version  of  Parseval’s  theorem 

”/ 


A^Sff  +  E-1]-1 

£  =  Sfv 


v=[  svv +  £;'], 


then 


Tr  [[St]"1]  =  Tr  [A]  +  Tr  A£  (r,  -  *'A0  1  £'A 


(3.29) 


The  second  tenn  can  be  simplified  by  applying  the  invariance  property  of  the  trace  under  cyclic 
permutation: 


A^-^AO^A 

=  Tr 

<>/-£' MY  '/aa£ 

_  JqxJq  JqxJq _ 

(3.30) 


=  Tr 


'(t/-^\Q-1(A0'A£ 

JqxJq  JqxJq  _ 


by  the  Hermitian  symmetry  of  A. 

(3.31) 


4  Summary  and  Plans 


In  this  report,  we  have  provided  a  brief  discussion  and  presentation  of  technical  details  of  the  work 
completed  in  the  previous  quarter  for  the  Video  Super-Resolution  Project.  General  plans  for  the 
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next  quarter  arc  the  following: 


(a)  Continue  development  of  two  mathematical/numerical  approaches  for  quantifying  information 
for  general  case  of  arbitrary  set  of  camera  perspectives  -  one  based  on  sparse  matrices  and  one 
based  on  aggregation  of  subimage  information  metrics. 

(b)  Carry  out  software  implementation  of  the  “aggregation"  information-theoretic  approach  (for 
computing  the  estimation  bounds  in  the  general  case)  in  (a)  and  do  preliminary  testing. 

(c)  Implement  and  test  software  implementation  of  the  methodology  for  adaptively  estimating  object 
content  for  the  information  metric  based  on  image  restoration. 

(d)  Significant  testing/validation  of  the  current  approaches  for  video  restoration  and  generation  of 
information  content  metrics  for  color/multispectral  image  data,  based  on  the  luminancc/chrominancc 
decompositional  approach  as  developed  in  the  last  quarter. 

(e)  Review  MPEG  noise  models  and  make  sure  that  some  realistic  version  is  included  as  part  of  the 
overall  information-theoretic  framework. 

(f)  Complete  write-up  of  full-perspective  approaches.  Do  a  staged-implementation  of  the  full  large- 
image  restoration  capability  w  ith  true  joint  estimation  of  the  w  arp  and  imagery,  and  in  particular 
the  joint  of  blur  and  perspective  implementation  -  carry  out  initial  implementation  of  this  joint 
object/warp  estimation  and  do  testing. 
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Appendices 


A  Mathematics:  Flat-earth  approximation  to  projective  trans¬ 
form 

For  discrete  pairs  of  images  with  significant  camera  rotation,  a  full  eight-parameter  projective 
transformation  is  required  to  align  the  images.  However,  for  closely  related  views  a  quadratic 
transformation  is  a  good  approximation  to  the  projective  transform  and  is  more  stable  to  com¬ 
pute!  Kumar).  Let  d  (.r,  y,  t x,t2)  be  the  displacement  field  describing  the  motion  between  images 
I  (.r,  y  t ! )  and  /  ( x ,  y,  t2).  Then,  the  image  at  time  tx  warped  into  alignment  with  the  image  at  time 
h  is 


The  quadratic  approximation  to  the  displacement  held  is  given  by 


dx  =  a  i  +  n2x  +  o3y  +  o7.r2  +  o8.t:y 
dy  =  o.i  +  o5.r  +  o6y  +  o7.ry  +  o8y2. 


The  implementation  of  the  quadratic  transform  involves  estimating  the  parameters  {o|,. . .  .n8}. 
creating  the  grid  defined  by  (dx,dy).  then  interpolating  the  image  at  time  f!  onto  that  grid.  For 
iterative  image  restoration  methods  it  will  likely  be  necessary  to  calculate  the  partial  derivatives  of 
the  quadratic  transform.  Since  the  quadratic  transform  is  an  interpolation  it  can  be  written  as  an 
operator 


I  ti,t2)  =  Q{t)I  {x,y,tx) 
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where  r  =  {ai, . . . ,  a8} 

0_ 

Oa7 


0  Q(t)I  =  ^  EE7  [7,'7,,]/'  (;r  (T)  -  »A.t,2/(t)  -  ™Ay) 


=  EE7[7'’ 


0 

Oa7 

Oh  O.r  Oh  Oy 
Ox  On  t  Oy  Oat 

Oh  Ox. 


=  ^2^/(7), //(]—/(  (x(t)  -  /iAj,  y  (r )  -  mAy) 

n  m 

=  EE7[71’ 


Ox  Oa7 


+EI>."i(s&)- 


Oy  Od, 


-  EE  ('i«.»i£)£+EE  ('!«.<)  &• 


So  the  partial  derivative  is  the  sum  of  two  images:  (1)  the  element-by-element  product  of  the 
original  image  interpolated  with  the  derivative  of  the  interpolation  kernel  and  the  derivative  of  the 


new  .r-gnd  with  respect  to  the  motion  parameter,  (2)  similar  to  (1)  but  replace  x  with  y . 


■7T-Q{t)  I 
(J(\  i 


EE 

tn  '  ' 


(T )  ,  NT^  V- 

da, 


(?/(T) 


mAu) 


Oy  ( r ) 

Ooj 


B  Parameterizing  Quadratically  Penalized  Spaces 

This  section  develops  a  natural  parameterization  for  quadratically  penalized  weighted  least-squares 
objective  functions,  and  uses  this  representation  to  develop  an  eflicient  means  of  searching  over  the 
regularization  parameter.  Consider  the  objective  function 


A  (y,  A)  =  (y-Bx)7  W  (y  -  Bx)  +  Ax7  Px  (B  1) 

with  the  solution  given  by 


x  —  argmin  A  (y.  A) 

.T 


(B.2) 


=  (BrWB  +  AP)  1  BTWy  (B.3) 

From  the  structure  of  the  solution  one  can  see  that  A  is  related  to  the  minimizer  x  through  the 
matrix  pencil  A  (Bi  WB.  P).  The  weighted  inner-product  can  be  eliminated  by  taking  its  Cholesky 
decomposition  and  defining  the  problem  over  a  rotated  basis 


Wh  :  W  =  W„7  W„ 
Bw  =  W|,B 


x=  (BhTBh  +  AP)",BIlTWlly 


(B.4) 
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As  a  result,  the  remainder  of  this  appendix  will  assume  W  =  I  without  loss  of  generality.  There 
exists  a  unique  rotation  of  the  space  such  that  (B7  B  +  AP)  1  is  congruent  to  a  diagonal  matrix 
which  is  a  function  of  A.  To  see  this,  consider  the  generalized  eigen-decomposition  of  B7  B  and  P 


3  V,  D  :  PV  =  BJ  BVD  (BV)5  BV  =  I 

(B.5) 

VT  (B;  B  +  AP)  V  =  I  +  AVrPV 

(B.6) 

=  1  +  AD 

(B.7) 

(BrB  +  AP)  1  =  VDAVr 

(B.8) 

Da  =  <Hag  ([1  +  ADn]  1 ,  •  •  •  ,  [1  +  ADPP]  ') 

(B.9) 

=  diag(dA) 

(B.10) 

from  which  the  estimate  may  be  restated  as 

x  =  VDA(BV)7  y 

(B.l  1) 

=  VDAa 

(B.12) 

=  V(dA  ©  a)  s.t.  a  =  ((y7  B)V)T 

(B.13) 

Now  consider  calculating  the  MSE  with  respect  to  this  parameterization 

MSE(x,  A)  =  ||x  —  x||2 

(B.14) 

—  cliaR  ( [  Du  (ADn  +  1)  , Dpp (ADPP  +  1)  J) 

(B.l  5) 

-TJ 

btO 

II 

(B. 16) 

[DAa]  =  dA  0  a 

(B. 17) 

(B.18) 

=  v^(D,a] 

(B. 19) 

=  v  (dA  0  a) 

(B.20) 

xTx  =  (DAa)7  VrV  (DAa) 

(B.21) 

3Da  bT*l  =  2VTV  (D»») 

(B.22) 

f}\  [X^]  =  (Da*)  V'V(dA0a)  +  (dA©a)  V'V  (DAa) 

(B.23) 

=  2(DAa)7  VTv(dA  3  a) 

(B  24) 

(B.25) 

=  (dA  ©  a)  Vx 

(B.26) 

Appendix  -  16 

Putting  this  all  together,  we  get  an  efficient  means  of  computing  both  the  MSE  and  its  derivative 
with  respect  to  A  in  terms  of  the  results  of  a  single  O  (n3)  calculation  to  find  the  canonical  space  (a 
generalized  eigen-value  decomposition). 


MSE  =  ||V(d\  ©a)  —  x||2 
^  [MSE]  =  -2  (VTV  (dA  0  a)  +  Vx) '  (dA  0  a) 


(B.27) 

(B.28) 
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1  Executive  Summary  /  Conclusion 

ASl  pursued  5  areas  of  study. 

Phenomenology  -  Modeling  and  measurements  were  analyzed  and  showed  that  HH 
polarization  was  preferred  and  gave  an  SNR  suitable  to  perform  the  desired  ISAR  image 
formation. 

System  modeling  -  System  modeling  was  done  to  establish  the  radar  requirements. 

Low  Power  Radar  RF  design  -  A  complete  strawman  transmitter  and  receiver  were  designed 
using  the  latest  generation  SiGe,  GaAs,  PHEMPT  technology  to  demonstrate  that  an 
extremely  small  radar  could  be  built  that  required  less  than  2  watts  to  operate  yet  meet  all 
performance  requirements. 

Low  Power  Radar  Digital  design  -  A  complete  strawman  radar  processor  was  designed 
using  current  generation  FPGA  technology,  cell-phone  RAM,  and  mobile  processors  to 
show  that  an  extremely  small  processor  could  be  built  that  operated  on  less  the  200mW. 

Antenna  design,  construction,  and  testing  -  An  antenna  small  enough  to  fit  the  intended 
application  was  fabricated  and  proved  to  have  all  the  performance  metrics  (bandwidth,  low 
dispersion,  bcamwidth,  and  small  size)  needed  to  perform  the  high  resolution  UWB  ISAR 
function  in  the  required  small  form  factor  system. 


TR- 1 1 090 1  Applied  Signals  Intelligence,  Inc.  FOR  OFFICIAL  USE  ONLY  (FOUO)  Page  6  of  49 


Subagreement  #100541Z1 


Riverine  Imaging 


September  30,  2011 


2  Phenomenology 

The  purpose  of  this  “quick  look”  analysis  was  to  understand  pro’s  and  con’s  of  vertical 
versus  horizontal  polarization,  and  why  one  might  be  better  than  the  other. 

2.1  Background  --  The  Model 

AS1  uses  an  internally  developed  a  fast  modeling  code  with  four  key  features.  It  models: 

•  Multipath,  including  1.  2,  and  3  bounce  terms 

•  Arbitrary  bandwidth/rcsolution  waveforms  (such  as  UWB  -  Ultra  Wide  Bandwidth) 

•  Both  SAR  (moving  radar)  and  ISAR  (moving  target)  cases,  and 

•  Arbitrary  ground/water  parameters  (e,  p,  a). 

Regarding  the  multipath,  it  includes  the  direct  (one  bounce)  path  (radar  to  target  and  back), 
2-bouncc  paths  (radar  to  target  to  ground  and  back  plus  radar  to  ground  to  target  and  back), 
and  3-bouncc  path  (radar  to  ground  to  target  to  ground  and  back).  Analysis  of  UWB  radar 
data  shows  that  these  paths  arc  by  far  the  most  dominant  and  accurately  model  UWB  radar 
data.  Regarding  the  SAR  and  ISAR  capability,  it  models  motion  along  an  arbitrary  path. 

2.2  Phenomenology  Summary' 

After  getting  model  results,  field  tests  were  made  that  confirmed  the  predictions  made  by 
compressed  UWB  mutipath  modeling — there  was  complete  agreement.  To  summarize: 

•  The  RCS  of  an  excessively  small  target  (small  empty  plastic  eanoc,  no  motor,  two 
people)  was  ~  0  dBsm  broadside  and  — 15  dBsm  at  end  of  a  SAR  aperture. 

•  The  radar  impulse  response  and  resolution  changes  significantly  over  the  range  swath 
due  to  the  significant  geometry  changes  in  the  one,  two,  and  three-bounee  multipath 
terms. 

•  IIH  is  superior  to  VV  due  to  combination  of: 

•  Clutter  from  water  ripple  is  suppressed  in  HH  polarization  due  to  multipath 
o  ~  25dB  @  5cm  height  and  >50dB  @  1  cm  height  at  100m  slant  range 

o  ~  20dB  @  5cm  and  45dB  @  1  em  at  50m  slant  range 

•  Target  (both  boat  and  contents  of  interest)  arc  enhanced  from  multipath  in  HH  pol. 

•  The  raw  target  RCS  is  larger  in  HH  -  Both  boat  and  contents  are  oriented  to  have  higher 
RCS  at  HH 

•  HH  experiences  consistent  enhancement  over  VV  for  targets  above  the  water 
o  With  lm  target  height 

■  HH  had  gain  except  from  58-78m  slant  range 

■  W'orst  case  loss  was  6  dB 
o  With  2m  target  height 

■  HH  gain/loss  curve  stretches  right  from  lm  ease 

■  HH  had  gain  except  from  108- 132m  slant  range 

■  Worst  ease  loss  was  6  dB 

2.3  Phenomenology  Details 
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Low  height  results,  representative  of  water  ripple,  are  shown  in  Figure  1,  while  lm  height 
and  2m  height  results  are  shown  in  Figure  2.  Figure  1  shows  there  is  a  strong  preference 
for  HH  polarization  to  reduce  the  clutter  caused  by  water  ripple.  Figure  2  shows  three 
important  features.  First,  that  the  typical  deep  nulls  that  would  happen  with  narrowband 
signals  is  mitigated  by  the  fact  that  the  UWB  waveform  resolves  multipath  at  the  shortest 
ranges  -  the  worst  ease  null  depth  is  reduced  to  only  6  dB.  Second,  that  there  is  signal 
enhancement  over  much  of  the  range-swath  of  interest.  Third,  that  this  outcome  is  caused 
by  that  fact  that  at  the  shortest  ranges,  the  multipath  is  resolved  by  the  UWB  waveform. 
At  other  ranges  the  multipath  is  in-phase  so  the  various  bounces  mostly  or  completely 
add  together. 
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Figure  3  plots  three  different  metrics  as  a  function  of  range  and  for  each  polarization:  (1)  the 
gain  over  the  direct  path,  (2)  the  received  pulse-width  (or  radar  resolution)  at  the  -lOdBe 
points,  and  (3)  the  peak  sidelobe  in  the  time-domain  (i.e.  down-range)  impulse  response. 
These  plots  highlight  and  provide  turther  insight  into  Figure  2.  Four  key  points,  labeled  A, 

B,  C,  and  D,  are  selected  for  closer  analysis.  Note  that  at  range  B,  there  is  a  convergence  of 
the  best  sidelobes,  a  peak  in  resolution,  and  a  peak  in  the  signal  enhancement  (or  gain  over 
the  direct  path). 
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Figure  3:  Correlation  between  gain,  resolution,  and  sidelobes  versus  range 
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Figure  4  shows  the  waveform  received  and  its  envelope  at  range  A,  B,  and  C.  Figure  5 
shows  the  received  resolution  “pulse  width”  and  sidelobe  levels  at  those  same  ranges.  At 
range  A,  the  1 -bounce  (closest  range),  2-bouncc  (middle  range),  and  3-bouncc  (farthest 
range)  terms  arc  all  resolved.  Note  that  the  2-bounce  return  is  tw  ice  the  amplitude  of  the 
other  terms.  This  doubling  is  because  there  are  two  paths  that  arc  identical  length  (i.c.  in- 
phase  at  all  frequencies)  that  add  together — radar-target-ground-radar,  and  radar-ground- 
target-radar.  At  range  B,  there  is  an  in-phase  convergence  that  gives  a  high  resolution,  high 
amplitude,  and  low  sidelobe  return.  At  range  C  the  return  is  just  the  reverse — low  amplitude, 
low  resolution,  and  high  sidclobcs.  Figure  6  shows  the  waveform  received  and  its  envelope 
at  range  D.  This  is  the  range  where  the  different  multipath  terms  add  in-phase.  Beyond  this 
point,  the  three  isolated  returns  seen  at  point-A  become  more  and  more  aligned,  and  the 
phase  of  the  center  term  becomes  1 80-dcgrccs  out  of  phase  with  the  other  terms,  causing  an 
ultra  wide  bandwidth  null.  To  summarize,  the  metrics  at  these  points  arc: 

•  A  -  40  dB  sidclobcs  and  multipath  resolved  (3  lobes) 

•  B  -  70  dB  sidclobcs  and  9  dB  gain  enhancement  due  to  multipath 

•  C  -  25  dB  sidclobcs  and  5  dB  loss  due  to  multipath 

•  D  -  Almost  12  dB  of  gain  -  (best  possible  with  4  paths  is  12  dB).  Note  that  it  is  12  dB 

as  opposed  to  6  dB  due  to  the  two-way  radar  path  as  opposed  to  a  one-way 
communications  path. 
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Figure  6:  Direct  &  Multipath  Received  Waveforms  at  Peak  at  R  -  134  m 

Measurements  carried  out  by  MTRI  showed  exactly  what  the  model  predicted.  (1)  There 
was  significantly  more  clutter  (water  ripple)  “noise”  in  VV  data  relative  to  HH  data.  The 
same  preference  for  HH  holds  with  trees  -  i.c.  tree  clutter  is  reduced  because  the  tree  RCS  is 
small  for  HH  and  large  for  VV.  (2)  HH  has  enhanced  target  (hidden  object)  RCS  because 
target  objects  arc  generally  oriented  along  hull  -  not  oriented  vertically.  In  other  words,  they 
have  long  horizontal  features  and  almost  no  vertical  features.  The  measurements  showed 
that  presence  of  metal  in  a  canoe  gave  higher  return  in  HH  data,  and  made  no  difference  in 
VV  data.  (3)  HH  has  enhanced  boat  RCS.  The  SNR  on  the  boat  was  significantly  higher  in 
HH.  The  boat  is  primarily  horizontal,  with  much  smaller  vertical  features.  Furthermore,  the 
water/hull  boundary,  both  far-side  and  near-side,  is  (a)  a  long  horizontal  feature,  and  (b)  a 
sharp  comer  at  the  water  line  where  the  hull  breaks  the  surface  of  the  water,  forming  a 
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depression  that  behaves  as  a  dihedral  “comer”  reflector,  with  an  especially  large  return  from 
the  far  side  of  the  boat.  The  watcr/hull  boundary  is  much  shorter  vertically.  Note  that  this 
watcr/hull  boundary  is  an  important  loud  return  since  it  is  very  predictable  and  can  be  used 
as  part  of  the  reference  system  for  an  ISAR.  The  HH  data  was  always  about  lOdB  louder 
than  VV  data.  The  worst  case  (small  empty  plastic  canoe,  no  motor,  two  people)  RCS  was  ~ 
0  dBsm  broadside  and  — 15  dBsm  at  end  of  a  SAR  aperture. 
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3  System  Modeling 

Objectives  of  this  effort  were  to  ( 1 )  derive  the  system  requirements  that  drive  the  radar 
design  such  as  minimum  and  maximum  PRF,  transmit  power,  noise  figure,  duty  cycles,  etc., 
and  (2)  estimating  the  power  budget  and  understand  radar  impaet  on  battery  life  by  (3)  doing 
a  strawman  design  to  support  the  power  and  size  budget  estimates. 

3.1  Overview 

Impulse,  poly-phase  eoded,  step-frequency,  and  ehirp  are  a  few  of  many  possible 
waveforms.  The  system  model  considered  the  time/power  versus  time  budget  for  three 
architectures: 

•  Impulse  -  Peak  power  &  volts  required  for  Impulse  -  PRF  limited  by  ambiguities 

•  Pulsed  Chirp  -  Chirp  length  limited  by  short  range,  PRF  limited  by  ambiguities 

•  Continuous  Chirp  -simultaneous  transmit  and  receive,  PRF  defined  by  sample  rate 
required  by  ISAR  and  target  motion  and  dynamics 

At  the  highest  level,  power  is  conserved  by  putting  blocks  to  sleep  (power  off).  Idle  modes 
(clock  off)  are  also  helpful,  but  even  when  clocks  arc  off,  static  power  will  kill  battery  life. 

In  other  words,  the  integration  of  a  small  power  over  a  long  time  still  comes  out  to  be  a  big 
number.  Integrating  zero  comes  out  zero  regardless  of  the  integration  time.  Therefore,  the 
key  to  low  power  operation  hinges  on  optimizing  sleep  time.  Because  there  is  usually  a 
significant  time  delay  between  a  power-off  state  to  a  fully  operational  state,  the  system  can’t 
use  many  short  sleep  time-slots,  but  must  use  a  few  long  sleep  time-slots. 

3.2  Background  on  chirp  and  stretch  processing 

Figure  7  illustrates  pulse  compression  theory  and  practice.  The  top  row  depicts  how- 
resolution  is  proportional  to  the  inverse  of  the  system  bandwidth.  In  other  words,  the  wider 
the  bandwidth  the  better  the  resolution. 

The  second  row  in  Figure  7  shows  how  pulse  compression  theory  works.  We  start  with  an 
impulse.  The  impulse  is  passed  through  an  all  pass  filter  H(co)  that  can  produce  a  chirp,  or 
some  phase  eoded  waveform.  This  waveform  gets  transmitted  and  received.  The  received 
waveform  is  passed  through  an  inverse  filter,  which  returns  the  signal  baek  to  the  impulse 
we  started  with.  The  purpose  for  doing  this  is  that  the  transmitted  waveform  ean  be  a  long 
duration  low  power  signal  with  low  peak-to-peak  volts.  After  compression,  it  ean  be 
equivalent  to  an  incredibly  large  peak-power  impulse  that  is  hard  to  build  due  to  the  high 
peak  voltages  and  currents  required. 

The  third  row  in  Figure  7  illustrates  the  case  where  H(co)  produces  a  chirp,  which  gets 
compress  to  an  impulse.  In  the  illustration,  the  chirp  has  a  1  GHz  bandwidth,  lasts  1  ms,  and 
gets  compressed  to  1  ns.  So  a  1  mW  chirp  is  compressed  to  the  equivalent  of  a  1  kW  peak 
power  transmitter. 
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A  key  feature  of  chirp  is  its  ability  to  use  what  is  called  “stretch”  or  “de-ramp”  processing. 
Figure  8  illustrates  stretch  processing.  The  top  half  of  Figure  8  is  a  plot  showing  the 
frequency  versus  time  of  the  transmitted  chirp  and  returns  from  targets  at  two  different 
ranges.  The  key  is  that  the  frequency  offset  between  the  transmitted  ehirp  and  the  returns,  is 
constant  and  at  a  different  frequency  for  each  range.  The  lower  part  of  Figure  8  shows  a 
receiver  structure  where  a  mixer  is  used  to  create  that  constant  frequency  that  is  a  function 
of  range.  All  the  offset  frequencies  are  captured  by  an  analog  to  digital  converter.  A  digital 
signal  processing  (DSP)  is  used  to  separate  the  various  frequencies  to  resolve  all  the 
different  target  ranges.  In  other  words,  the  DSP  compresses  the  signal  back  to  what  would 
have  been  received  had  an  impulse  been  transmitted  and  received  instead  of  the  ehirp.  The 
DSP  requires  a  Fourier  transform  and  a  de-skew  operation  to  make  the  output  identical  to  an 
impulse  based  system. 

Why  is  stretch  important?  Because,  when  the  compression  ratio  is  high,  the  ADC  bandwidth 
does  not  need  to  eover  the  entire  bandwidth  of  the  transmitted  signal.  Instead,  the  ADC 
bandwidth  is  proportional  to  the  “range  swath”  or  the  distance  between  the  shortest  range 
and  farthest  range  of  interest.  More  precisely: 


ADC  Sample  Rate  (Bandwidth)  oc  (Max  Range  -  Min_Range)*Ramp_Ratc  (MHz/ps)  or 

ADC  Sample  Rate  (Bandwidth)  oc  (Max  Range  -  Min_Range)*Bandwidth/(Chirp  Duration) 

Since  ADC  power  goes  up  with  sampling  frequency,  maximizing  the  chirp  duration 
minimizes  the  ADC  power.  Continuous  chirp,  i.e.  transmitting  and  receiving  simultaneously, 
maximizes  the  chirp  duration.  It  minimizes  the  average  power  drawn  while  the  radar  is 
active.  Where  docs  stretch  break  down?  It  breaks  down  if  the  compression  ratio  is  too  small. 
Indeed,  it  ean  even  make  the  ADC  bandwidth  larger  than  the  signal  bandwidth.  Due  to  end 
effeets,  it  is  most  efficient  when  the  compression  ratio  is  high  -  or  in  other  words,  where  the 
chirp  is  very  long  relative  to  the  compressed  pulse. 
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Pulsed  chirp  allows  the  receiver  to  operate  while  the  transmitter  is  off.  A  pulsed  chirp 
system  generally  has  a  huge  advantage  in  dynamic  range  since  the  receiver  is  not  trying  to 
operate  while  the  transmitter  is  potentially  blinding  the  receiver.  In  long  range  applications, 
it  is  the  only  way  to  go.  But  it  drastically  raises  the  battery  drain  in  short  range  applications. 
In  this  ease,  because  of  the  short  range,  the  chirp  cannot  be  long.  Therefore,  the  efficiencies 
of  high  compression  ratio  ehirp  disappear.  Even  if  ehirp  portions  were  transmitted  so  that  an 
entire  chirp  is  sent  at  about  the  same  rate,  the  chirp  rates  and  ADC  clock  go  up  by  more  than 
a  faetor  of  2.  The  air  between  the  radar  and  nearest  range  bin  must  be  filled  with  ehirp 
energy,  with  enough  time  for  the  transmitter  to  be  turned  off,  the  receiver  turned  on,  and  all 
the  switching  transients  to  settle,  so  that  the  receiver  can  be  receiving  the  energy  from  the 
nearest  range.  The  receiver  must  then  stay  aetive  long  enough  to  receive  the  energy  to  the 
farthest  bin.  This  leads  to  the  >2X  faetor.  Beyond  this  core  faetor,  the  hardware  required  to 
support  generating  tightly  controlled  phase  on  the  chirp  sections,  with  fast  settling  times, 
adds  complexity  and  draws  all  the  more  power. 


Tx  Chirp  ...  ,  _  .  _ 

Near-far  Dynamic  Range  >  120  dB 


Figure  8:  De-Ramp  (“Stretch”)  Implementation 

The  combination  of(l)  the  reduced  dynamie  range  requirement  of  short  range,  (2)  the  short 
range  swath,  and  (3)  the  high  dynamie  range  obtainable  with  today's  semiconductor 
technology,  i.e.  mixers,  amplifiers,  and  ADC,  make  the  use  of  continuous  chirp  a  viable  and 
lower  power,  option. 

3.3  Geometry  and  sampling  (PRF  and  number  of  points)  calculation 

Figure  9  shows  a  side  view  of  the  radar/target  geometry.  Figure  10  shows  a  top  view  of  the 
geometry  along  with  a  spreadsheet  that  was  used  to  look  at  the  PRF  requirements,  the 
duration  of  an  aperture,  and  the  number  of  pulses  required.  The  final  highlighted  columns 
show  the  number  of  samples  in  an  aperture  assuming  a  250  Hz  PRF,  which  gave  some 
margin  over  the  208  Hz  required  for  a  worst  ease  wide  integration  angle  at  close  range.  It 
reduces  smearing  caused  by  the  range  sampling  widow  occurring  over  a  4ms  time  span. 
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R  =  22  -  200  m,  0  =  3°  -  26°  @  20m  standoff  from  bank 
R  =  10.2  -  200  m,  0  =  3°  -  78°  @  2m  standoff  from  bank 


Figure  9:  Scenario  Geometry 
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Figure  10:  Energy  Calculation 


3.4  Transmitter  Power  Calculations 

Range  compression  gain,  ISAR  aperture  gain,  and  range  and  cross-range  resolution  are  a 
function  of  the  window  taper  used  on  each.  Figure  1 1  lists  the  loss  factors  for  various 
window-tappers  assuming  Nyquist  sampling  -  no  up-sampling. 
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Figure  1 1:  Characteristics  of  various  window  taper  functions 


Figure  12  shows  a  calculation  of  system  losses.  The  loss  for  detections  uses  only  the  first 
two  rows  since  detection  does  not  involve  imaging. 


Loss 

(dB) 

Description 

1.4 

Range  Compression  Kaiser  Window  0=5.44  (~samc  as  Hamming) 

1.7 

Scalloping  loss  (reduces  to  .41. 1  dB  with  2X/4X  zero-pad 
oversampling) 

1.4 

Kaiser  for  aperture  (imaging) 

1.7 

Scalloping  loss  for  Imaging 

6.2 

Total 

3.1 

For  Detection 

Figure  12:  System  losses 


Figure  13  defines  the  parameters  and  lists  the  equations  used  for  the  transmitter  power 
calculation. 
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Figure  13:  Equations  used  to  compute  transmit  power  requirements 


To  conserve  power,  the  system  runs  in  two  modes.  A  detection  mode  with  a  high  sleep  time 
duty  cycle.  In  this  mode,  the  system  wakes  up,  transmits  a  single  pulse,  and  determines  if  a 
target  is  present  that  should  be  imaged.  If  not,  the  system  goes  back  to  sleep.  If  a  target  is 
detected,  the  radar  goes  into  the  second  mode,  where  it  sends  enough  pulses  to  image  the 
target. 

Figure  14  lists,  for  imaging  an  interior  object,  the  parameter  values  that  went  into  the 
equations  in  Figure  13. 
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Figure  14:  Parameter  values  used  for  transmit  power  calculation 


Using  these  parameters,  the  transmitter  requirements  are: 
Et  =  0.34  mJ  @  1  GHz 
=  0.68  mJ  @  1.5  GHz 
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Assuming  the  same  power  is  used  for  single  pulse  detection  in  the  detection  mode,  i.e.  n  =  1 
(single  pulse  detection),  Lp  =  1  (no  penetration),  and  Ls  =  .5  (3  dB  system  loss)  we  find  o  - 
1.6  dBsm 

3.5  Case-1  -  Tx  &  Rx  in  separate  time  slots 

Figure  15  lists  the  calculations  used  to  look  at  a  short  chirp  and  impulse  based  transmitter. 
Note  that  the  impulse  system  power  is  approximately  the  same  voltage  as  used  on  the  Army 
Research  Lab  (ARL)  BoomS AR.  At  this  power  level,  the  pulse  generator  is  large  and 
inefficient.  While  a  higher  PRF  might  be  used  to  reduce  the  voltage  requirements,  doing  so 
has  two  drawbacks.  First,  a  higher  PRF  also  requires  a  higher  duty  cycle  on  the  ADC  and  the 
processing  that  must  digest  the  increased  data  flow.  Second,  a  higher  PRF  ultimately  runs 
into  an  ambiguity  issue.  The  ADC,  for  both  the  50ns  and  2  ns  pulse  needs  to  cover  the  entire 
transmitted  bandwidth.  Currently,  8-bit  1.5  Gsps  ADCs  burn  6  Watts  and  an  8-bit  2.2  Gsps 
ADCs  burn  8  Watts.  In  contrast,  the  16-bit  20  Msps  ADC  selected  for  the  continuous  system 
burns  only  60  mW. 
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p-p  into  50  ohms  @  250  Hz 

Figure  15:  Peak  transmit  power  eale,  Tx  &  Rx  in  separate  time  slots 


3.6  Case-2  0  Tx  &  Rx  simultaneously 

In  this  ease,  we  minimum  the  peak-power  and  the  ADC  bandwidth  by  operating  at  a  100% 
duty  cycle.  Given  the  250  Hz  PRF,  we  can  use  a  4ms  chirp  -  (125  Hz  Triangle  Wave,  4ms 
up-chirp,  4ms  down-chirp).  The  transmit  amplifier  only  needs  to  make  0.68mJ/4ms  = 

1 70mW  or  +22dBm  (7.9  V  p-p  into  50  ohms).  This  low  power  level  has  two  major  impacts. 
First,  small  geometry  wide  bandwidth  and  efficient  transistors  in  GaAs,  InP,  GaN,  and  SiGc 
can  easily  withstand  this  voltage  level.  Efficiency  at  higher  power  levels  typically  depends 
on  narrowband  tuning.  A  wideband  system  does  not  have  this  luxury,  and  is  therefore 
usually  much  less  efficient.  The  viability  of  using  small  geometry  transistors  will  raise  the 
efficiency.  Second,  the  receiver  has  enough  dynamic  range  to  manage  this  relatively  low 
transmitter  power  level  -  allowing  the  receiver  to  operate  while  transmitting. 

Generation  of  this  4ms  long  chirp  is  significantly  easier  (smaller,  lower  power,  less 
expensive)  than  a  50ns  chirp.  Furthermore,  it  is  easier  to  control  so  as  to  extend  the  battery 
life.  For  example,  the  power  can  be  reduced  depending  on  either  (a)  range  -  i.e.  w  hether  the 
target  is  detected  at  the  far  side,  or  close  side  of  the  river,  or  (b)  the  detected  SNR,  in 
combination  with  the  detected  speed  of  the  target.  The  4  ms  time  is  long  enough  that  the 
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pulse  length  could  be  reduced.  For  example,  the  pulse  could  be  reduced  to  2  ms  and  the 
system  could  sleep  for  2ms.  This  would  save  the  battery  yet  not  change  the  sample  rate 
(PRF)  that  the  1SAR  needs.  If  the  target  is  moving  slower,  then  the  radar  can  reduce  the 
power  by  reducing  the  PRF  and  sleeping  between  pulses.  With  this  approach,  the  amplifiers 
remain  biased  at  their  most  efficient  operating  point.  Another  power  reduction  option  is  to 
bypass  and  switch  off  the  final  TX  power  amplifier. 
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4  RF  Hardware  Strawman  and  Power  Evaluation 

4.1  Exciter  Strawman  Hardware 

The  exciter  needs  to  make  an  extremely  linear  chirp  in  order  to  achieve  a  reliable 
compression  ratio  of  2*  10°.  Common  approaches  are: 

•  Linearize  a  voltage  controlled  oscillator  (VCO)  by  driving  it  with  an  ADC  via  a  look-up 
table  (LUT). 

This  open-loop  approach  is  not  viable  as  it  cannot  achieve  the  linearity  required  to 
aehieve  the  needed  compression  ratio. 

•  Direetly  make  the  chirp  with  a  DDS  (direct  digital  synthesizer)  -  a  digitally  produced 
waveform  driving  an  ADC. 

This  approach  has  the  advantage  of  nearly  instant  wake-up/sleep  and  no  VCO  drift, 
warm  up,  settling  time,  or  loop  stability  issues.  It  is  not  viable,  however,  due  to  power 
required  to  get  the  required  clock  rate  -  on  the  order  of  a  few  Watts. 

•  Make  low  bandwidth  chirp  with  a  DDS  and  lock  a  VCO  using  a  standard  PLL  to 
multiply  the  low-BW  chirp. 

Figure  17  illustrates  a  strawman  design  using  the  latest  generation  low  power  parts  to 
implement  this  approach.  Advantages  are  that  the  power  is  only  ~  0.63  Watts  and  the 
PLL  loop  bandw  idth  is  wide  so  tracking  at  endpoints  is  fast  and  VCO  linearization  is 
finely  spaced.  Disadvantage  is  that  the  phase  noise  of  the  DDS  is  multiplied.  As  an 
example,  with  M=16  the  DDS  might  generate  a  50MHz  to  100  MHz  chirp,  and  the 
VCO  would  generate  an  800  MHz  to  1 600  MHz  chirp. 

•  Make  low  bandwidth  chirp  with  a  DDS  and  lock  a  VCO  using  a  fraetional-N  PLL 
eonneetion. 

Figure  16  illustrates  a  strawman  design  using  the  latest  generation  low  power  parts  to 
implement  this  approach.  Advantages  are  that  the  power  is  only  ~  0.63  Watts,  the  DDS 
output  is  filtered  around  the  reference  frequency  providing  some  noise  reduction,  the 
DDS  runs  slower  and  those  slower  running  portions  of  the  chip  will  draw  less  power. 
The  disadvantage  is  that  the  settling  time  may  not  be  fast  enough. 

•  Make  low  bandwidth  chirp  with  a  DDS  and  direetly  multiply  it  up  to  the  required 
frequency  range.  This  has  the  advantage  of  allowing  the  fastest  power  cycling  and 
extremely  fastest  settling  times.  The  multiplication  factor  could  be  low  so  that  a  single 
amplifier  eould  be  used  to  boost  the  desired  harmonic.  Power  used  in  the  VCO  and 
amplifiers  in  Figures  16  and  17  would  be  used  for  the  RF  amplifier.  The  disadvantage  is 
that  high  fundamental  rejection  (i.e.  good  balance)  is  needed  in  the  multiplier  stage. 

The  good  news  is  that  the  last  three  approaches,  and  in  particular,  the  diagrams  in  Figures  16 
and  1 7,  show  that  the  required  low  power  level  is  achievable. 
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Figure  16:  VCO  linearized  &  locked  in  a  fraetional-N  (DDS)  configuration 


Ref  OSC 

400MHz 
150ppm 
-20°to  80°C 
3V@1.2mA 
100  mW 

2-1 0ms 
Startup 
2x  .13  in 
CVT32-10 
CSX532T 


Clock  Gen  Chip 


AD9540 

400 mW 


£ 


tM  r>  pfd 

Phase 
Freq 
Det 


DDS  >_ 
fN-^DAC  ) 
4A(t) 


LPF 

OpAmp 

25mW 


VCO 

1 

.72-1.75  GHz 

1 

1 

.99-1.815  GHz 

1 

1-2  GHz 

1 

1 

5V@20mA  , 

1 

100  mW 

1 

V585ME41LF 

1 

5dBm  out 

1 

1 

,5x  .5  in 

1 

1 

Amp 

15  dB  gain 
140mWout 
21  5dBmcp 
5V@60ma' 
300  mW 
PMA-5456 


Y 

Filt- 


Cpl 


Amp 

alETH1 


1 


To 

Rec 


10  dB  gain 
IWout 
30dBm  cp 
5V@500ma 

2.5W 

hmc45lst89 


Figure  17:  VCO  Linearized  &  Locked 


4.2  Receiver  Strawman  Hardware 


Figure  18,  shows  a  strawman  design  for  the  receiver.  It  uses  a  classic  direct  conversion 
approach.  This  approach  eliminates  the  second  chirp  generator  that  is  often  used  to  offset  the 
range  swath  and  center  the  baseband  signal  to  be  at  DC  in  the  center  of  the  swath.  Doing  so 
saves  all  the  power  of  another  Figure  16  or  Figure  1 7  block  except  for  the  crystal  oscillator. 
High  dynamic  range  is  achieved  by  using  the  latest  generation  PHEMPT  based  RF 
amplifiers,  triple  balanced  mixers,  SiGe  baseband  amplifiers,  and  16-bit  ADCs.  Very  high 
resolution  (dynamic  range)  and  high  speed  ADC’s  are  available  at  low  cost  and  low  power, 
as  shown  in  Figure  18.  The  relatively  high  20  MHz  sample  rate  of  the  converter  shown,  can 
be  used  to  reduce  the  duty  cycle.  Any  final  design  will  need  to  consider  the  power  tradeoffs 
between  the  ADC  speed/power  and  the  power  savings  possible  via  a  reduction  in  the  duty 
cycle. 
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Figure  16,  Figure  17,  and  Figure  18  take  advantage  of  the  latest  generation  parts  to  achieve 
low  power.  For  example,  the  ultra  low  power  crystal  oscillator  developed  was  developed  for 
cell  phone  &  high  volume  portable  applications.  The  high  power  output  amplifier  is 
configured  with  GaAs  switches  to  allow  bypassing  an  allow  adaptive  changes  in  the  duty 
cycle  to  enhance  battery  life.  We  believe  these  diagrams  that  use  real  parts  give  strong 
support  for  battery  life  numbers  and  radar  performance  numbers. 
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Figure  18:  Receiver 


4.3  RF  Hardware  Power 


At  100%  duty  cycle,  the  transmitter/exeitcr  draws  930mW  and  the  receiver  draws  630  mW, 
for  a  total  of  1 .56W.  Because  of  bias  currents  often  dominating  power  drain,  typically 
adjusting  output  power  often  barely  affects  battery  life.  To  provide  a  path  to  extend  the 
battery  life,  we  will  adjust  the  pulse  durations  via  duty  cycle  to  reduce  both  transmit  and 
battery  power  simultaneously.  With  a  sample  rate  of  20  MHz,  and  assuming  1 ,5X  Nyquist 
sampling,  the  duty  cycle  can  go  from  1 00%  (4ms)  to  1 .5%  (0.06ms).  This  is  based  on  the 
frequency  that  the  ADC  must  capture  is  2  K  Rmax  /  C,  where  K  is  the  Hz/Sec  ramp  rate  of  the 
ehirp,  and  C  is  the  speed  of  light.  Assuming  1 ,5x  Nyquist  Sampling,  the  ADC  sample  rate  Fs 
=  3  Fo  (f=  6  K  Rmax  /  C,  and  the  maximum  ramp  rate  Kmax  =  Fs  C  /  (6  Rmax)-  So  at  20  Msps, 
the  chirp  rate  can  be  as  high  as  600MHz/0.06ms. 

Assuming  R4  power  sealing,  targets  at  34m  or  less  can  use  1 .5%  duty  cycle.  The  duty  cycle 
=  (7?/100m)4.  Figure  19  plots  the  battery  watts  consumed  versus  target  range,  during 
imaging.  It  assumes  a  minimum  RCS  target.  A  higher  RCS  target  could  use  less  battery  - 
the  entire  curve  moves  down.  The  knee  in  the  curve  is  due  to  20  Msps  sample  rate  driving 
the  minimum  duty  cycle  and  the  high  PRF  needed  to  capture  the  higher  ISAR  angular  spin- 
rate  at  short  ranges. 
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Slant  Range  m 

Figure  19:  Battery  Power  During  Imaging 
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5  Processor  Strawman  Design  and  Power  Evaluation 

5.1  Overview 

The  purpose  of  this  task  was  to  derive  accurate  power  budget  so  there  was  a  better 
understanding  on  how  the  complete  radar,  both  RF  and  Digital  sections,  impacts  battery  life. 
A  key  part  of  this  work  was  the  derivation  of  a  set  of  modes  (like  sleep,  blink,  run,  detect, 
image)  and  a  system  time  line  or  operational  concept  that  sequenced  between  these  modes 
so  as  to  minimize  power.  The  derived  power  numbers  are  supported  by  mapping  the 
required  processing  onto  real  hardware  -  in  this  case,  a  current  generation  Altera  FPGA, 
cell-phone  RAM,  and  a  cell-phone  processor,  etc.  The  processing  captured  the  key  signal 
processing  steps: 

•  Range  compression 

•  CFAR  detector 

•  (ok  (range  migration)  processing 

•  Feature-creation  (length,  w  idth,  RCS,  etc.) 

Ignored  processing  included  scaling/packing/comprcssing  of  the  data  for  transmission,  and 
trajectory  estimation.  While  the  trajectory  estimation  algorithm  was  ignored,  the  processing 
load  will  likely  be  smaller  after  its  insertion.  This  reduction  is  anticipated  because  the 
estimated  trajectory  will  allow  decimation  of  data  set  prior  to  FFT  processing. 

After  a  survey  of  DSP  and  FPGA  options,  we  chose  an  FPGA  approach  for  baseline,  and 
used  an  Altera  Cyclone  III.  Why  the  shift  away  from  CPU’s  and  toward  FPGAs  for  DSP? 
Until  recently,  IC  scaling  has  been  on  both  density  and  speed.  But  today,  as  the  feature  sizes 
shrink  further,  the  IC  performance  scales  primarily  in  density,  and  NOT  speed.  This  is  a 
major  shift/driver  on  semiconductor  industry.  An  example  outcome  is  the  emergence  of 
multi-core  processors  as  the  “standard”  processor. 

DSP  chips  do  have  their  advantages.  The  “C”  programmability  of  DSP  is  a  huge  advantage 
for  DSP  chips.  It  makes  the  development  cycle  is  also  much  more  manageable  since  there  is 
a  large  base  of  engineers  who  know  “C”,  it  is  much  more  easily  documented  and 
maintained,  and  it  allows  fast  compiles  so  that  test-fix-test  iteration  loops  close  faster.  There 
arc  no  “routing  and  timing  closure”  issues.  For  all  these  reasons,  DSP  chips  arc 
advantageous  when  control  is  a  major  part  of  its  job.  FPGAs,  however,  take  better  advantage 
of  the  density  sealing  and  can  do  more  processing  in  a  smaller  size  and  with  lower  power. 
Figure  20  lists  why  this  outcome  is  the  case. 


TR- 1 1 090 1  Applied  Signals  Intelligence,  Ine.  FOR  OFFICIAL  USE  ONLY  (FOUO)  Page  25  of  49 


Subagreement  #100541Z1 


Riverine  Imaging 


September  30,  2011 


DSP 

FPGA 

Pipeline  —  Fixed  1 -pipeline  fits  all 

Requires  blazing  speed  on  short  fixed 
pipeline 

Can’t  change  balance  between  different 
functions 

(c.g.  DSP  chip  has  predefined  number  of 
adders  and  multipliers) 

Pipeline  --  Custom  pipelining  on  everything 

Takes  advantage  of  today’s  sealing  of 
density 

Capacity  to  have  many  parallel  functional 
blocks 

Can  choose  more  adders  and  fewer 
multipliers  or  vice-versa 

Highly  serial,  each  functional  block  is 
optimized  to  the  transistor  &  blazing  fast 

Highly  parallel,  functional  blocks  can  be 
slower  because  there  arc  so  many  more  of 
them 

Multiplexed  Bus  Centric 

Bus  shared  by  multiple  blocks 

Lots  of  connections-lots  of  stray  C 
Multiplexing  causes  toggling  that  wastes 
power 

Point-to-Point  Bus  Centric 

Most  ‘‘busses”  are  point-to-point 

Low  number  of  connections  -  low  stray  C 

No  multiplexing,  if  data  docs  not  chance, 
no  power  is  wasted  charging/diseharging  C 

1  clock  fits  all 

No  granular  control  of  what  is  getting  a 
eloek 

no  granular  control  of  what  runs  fast  and 
what  runs  slow 

Multiple  Clock  domains  easy 

Fine  control  over  what  is  clocked  &  when 
Takes  fine  advantage  of  duty  cycle  on 
every  block 

Fine  control  over  what  runs  fast  and  slow 
Examples: 

Memory  controller  rate  governed  by  memory 
Function  blocks  run  at  rates  limited  by  their 
own  logic 

Interface  logic  runs  at  rate  governed  by  the 
interface 

General  purpose  nature  results  in  powered 
and  eloeked  blocks  that  arc  not  used  for  our 
design 

Less  “stray  C”  because  unused  logic  is  not 
eloeked/eonneeted 

If  a  section  cannot  be  powered  down, 
it  burns  static  power  but  is  not  eloeked 

Already  hardened  -  very  optimized  chip 
layout 

Power  &  Performance  can  be  enhanced  by 
hardening  the  design  -  e.g.  “HardCopyTM” 
where  FPGA  code  is  translated  into  an 
optimized  chip  (no  longer  programmable)  - 
like  a  custom  pipelined  DSP. 

Figure  20:  DSP  versus  FPGA  drivers 

After  a  survey  of  memory  options,  we  ehose  eell  phone  memory  for  the  baseline.  The 
combination  of  low  power  and  high  speed  for  this  memory  is  a  major  breakthrough  for 
enabling  this  application. 

After  a  survey  of  ultra  low  power  processors  for  housekeeping  functions  and  the  low  power 
deteet-modc  processing,  we  ehose  a  “mobile”  ARM,  again  designed  for  low-power  portable 
consumer  electronics  products. 
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We  designed  a  pipelined  FPGA  architecture  to  implement  the  required  processing.  We 
implemented  with  Altcra’s  design  tool,  the  most  power  consuming  (many  gates)  functions 
(e.g.  an  FFT  block,  etc.)  to  get  reasonably  accurate  size  and  power. 

The  system  has  two  modes  -  detection  and  imaging.  Each  has  separate  computation  path. 

Detection  runs  at  a  1  Hz  rate  to  achieve  ultra  low  power.  Once  per  second,  the  ARM 
processor  wakes  up,  powers  the  radar  up  long  enough  to  collect  one  or  a  few  pulses, 
processes  the  full  100m  range  of  data  and  cither  (a)  updates  the  background  if  no  target  is 
detected,  or  (b)  wakes  up  the  FPGA  and  begins  making  an  ISAR  image.  This  mode  is 
designed  to  minimize  power  by  maximizing  the  sleep  time. 

Imaging  runs  at  250  Hz  radar  pulse  rate.  To  save  power,  this  mode  only  runs  when  a  target 
is  present.  The  background  is  not  updated  since  targets  arc  in  view.  This  mode  is  designed  to 
minimize  the  powered-up  time  of  the  FPGAs.  The  processing  is  split  between  two  FPGAs  so 
that  when  back-to-back  images  are  needed  one  can  still  have  a  sleep  cycle,  resulting  in  lower 
overall  power.  The  output  product  can  take  several  forms  suitable  for  enhanced  processing. 
These  include  a  range  compressed  data  set,  a  simple  focused  data  set  that  is  invertible,  an 
auto  focused  data  set  that  is  not  invertible,  and  features  like  length,  width,  RCS,  and  relative 
locations  of  largest  peaks  and  their  RCS. 

5.2  Digital  System  Hardware  Design 

Figure  21  shows  a  block  diagram  of  the  digital  hardware  that  includes  the  processing 
system.  Figure  22  lists  the  features  and  functions  of  the  selected  “mom”  processor.  Figure 
23  lists  the  features  and  functions  of  the  selected  memory.  Figure  24  lists  the  features  and 
functions  of  the  selected  FPGAs. 


Figure  21:  Block  Diagram 
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ARM  Cortex  M3  -  (e.g.  STM32F103) 

Functions 

Mobile  Oriented 

Mom  functions  -  good  general  epu 

35ma@2.7V  (lOOmw) 

Entire  Detection  Mode 

70ma  with  all  peripherals 

—  Processing  (Stretch,  CFAR) 

5uA  deep  sleep 

—  Background  memory  processing 

watch-dog  timer  ( 1  see  wake  up  &  detect) 

Flash  Memory  management 

Single  cycle  32-bit  multiply 

Power  management 

2  cycle  MAC  (multiply  &  accumulate) 

turning  on/off  FPGA,  radar,  maybe  eomm’s. 

64kB  on-board  memory 

Figure  22:  ARM  Cortex  M3  “Mom”  processor 


Cellular  RAM  -  MT45W8MW16  =  8Mxl6 

Functions 

Also  called  PSRAM  (Pseudo  Statie  RAM) 

Memory  to  collect  full  2D  aperture  data 

<10  uA  power  down 

Memory  for  2D  FFT  during  processing 

200uA  standby 

Combined  Write  &  Read  on  8k  by  lk  256ms 

35ma(a),l  .8V  =  full  speed 

Serateh  pad  and  Background  Reeord 

10/35  ns  (in-page/other  page) 

Figure  23:  Memory 


FPGA  -  Altera  Cyclonelll  LS 

Functions 

F.P3CLS70 

FPGA-1  -  Imagine  mode  only  (off  in  detect) 

Has  200  18x18  multipliers 

ID  processing  up  to  first  exfiltration  point 

8k  FFT  in  8k  eloeks 

FPGA-2  -  Imagine  mode  only  (off  in  detect) 

~384  kB  memory 

2D  processing  beyond  first  exfiltration  point  i.e.  past  first 
“corner  turn” 

60ms  power  up  &  eonfig  time 

Powered  off  between  output  of  2D  target  image  chips 

150mW  statie  power 

Up  to  8  targets  (e.g.  different  ranges)  can  be  processed  by 
reducing  sleep  time 

Figure  24:  FPGAs 


While  the  two  FPGAs  could  be  combined  into  1  Larger  FPGA,  the  split  allows  FPGA-2  to 
be  shut  off  between  images.  It  is  fast  enough  to  finish  an  image  before  a  new  aperture  has 
been  eolleeted.  This  faet  allows  two  benefits.  First  statie  power  is  conserved  when  there  are 
baek-to-baek  targets.  Second,  the  “excess  capacity”  can  be  used  to  form  image  chips  of 
multiple  targets,  such  as  two  targets  coming  down  the  river  together.  They  might  share  the 
same  aperture,  but  be  at  different  ranges. 

5.3  Main  Power  Consuming  DSP  Blocks 

The  FFT’s  used  block  floating  point  to  achieve  a  high  dynamic  range  fixed  point 
implementation.  There  are  three  classes  of  implementations  for  fixed  point  add/multiply. 

•  (1)  Shift  input  bits  so  that  it  is  impossible  to  overflow.  For  example,  use  only  9  bits  out 
of  16  for  a  128  point  FFT.  Clearly  this  choice  results  in  a  huge  loss  of  dynamic  range. 

•  (2)  Fixed  Sealing  -  Perform  a  2-bit  shift  after  each  butterfly  stage. 

•  (3)  Block  Floating  Point  -  Perform  a  bit  shift  only  if  needed  after  each  butterfly  stage. 
This  method  requires  output  of  an  exponent  to  know'  how  data  was  sealed. 
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Figure  25  illustrates  the  difference  in  the  noise  floor  between  fixed  scaling  and  block 
floating  point.  Block  floating  point  is  >20dB  better. 


Time  -  ps  time  -  ps 

Block  Floating  Point  pulse  Fixed  scaling  pulse 

Figure  25:  Reason  For  Preference  For  Block  Floating  Point  FFT 


The  power  budget  is  dominated  by  large  blocks  -  i.c.  FFT  blocks.  The  term  LE  is  Altera’ s 
term  for  “Logic  Element”  which  represents  a  block  that  can  do  any  4-input  combinatorial 
logic  with  or  without  a  register.  One  large  block  is  an  8k  x  lk  2D  FFT.  This  block  has  the 
following  features: 

•  Allows  variable  size  (max  8k) 

•  Uses  -9000  LEs 

•  Uses  96  1 8-by- 1 8  Multipliers 

•  Uses  64  block  RAMs 

•  Pipelined  so  that  there  is  1  complex-sample  in  &  out  per  clock 

•  ~300mW@  100MHz 

•  2  *  8k  *  1  k  cycles  =  1 60ms  with  fast  RAM 

•  It  takes  ~256ms  to  complete  because  cell  phone  RAM  acecss  time  dominates 
Higher  speed  memory  is  not  viable  due  to  the  power  dissipation. 

Another  large  block  is  the  4k  FFT  that  must  be  done  1600  times.  This  block  has  the 
following  features: 

•  Uses  -6000  LEs 

•  Uses  48  1 8-by- 18  Multipliers 

•  Uses  76  block  RAMs 

•  Pipelined  so  that  there  is  1  complex-sample  in  &  out  per  clock 

•  -200mW@  100MHz 

Miscellaneous  other  functions  are  small  enough  to  neglect  relative  to  the  huge  6000  LE  and 
9000  LE  FFT  blocks.  These  include: 

•  Square  root:  -400  LEs,  pipelined  -  1  answer  per  clock  -  only  used  on  1  vector  per  image 
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•  Interpolate:  pipelined,  similar  to  N-point  FIR  1  in  &  1  out  per  clock 

•  Where  N  is  the  number  of  input  samples  covered  by  resampling  filter 

•  This  is  the  largest  block  behind  the  FFTs 

•  Zero  pad:  1  clock  pcrsample  (zero  padding  is  always  to  internal  memory) 

•  Vector  multiply:  pipelined  1  in  &1  out  per  clock  regardless  of  real  (1  mult)  or  complex 
(4  mult) 

•  2D  Hanning  Window  multiply  -  1  clock  per  point 

•  vectorized  so  it  is  always  running  in  1 D  in  internal  memory 

•  CFAR:  pipelined,  1  in&l  out  per  clock 

5.4  CFAR 

Figure  26  shows  a  standard  cell-average  constant  false  alarm  rate  (CA-CFAR)  detector.  This 
structure  can  be  pipelined.  An  example  FPGA  pipelined  design  is  in  “Design  and 
Implementation  of  a  CFAR  Processor  for  Target  Detection”  by  Cesar  Torrcs-Huitzil,  Rene 
Cumplido-Parra,  Santos  Lopez-Estrada,  Computer  Science  Department,  INAOE,  Apdo. 
Postal  51  &  216,  Tonantzintla,  Puebla,  Mexico.  The  summary  listed  below  show  how  the 
block  is  tiny  relative  to  the  FFTs. 

•  Process  is  pipelined  to  I  in  &  1  output  per  clock 

•  Bit  w  idth  depends  on  using  approximate  Mag  function  or  Mag'  function 

•  Assume  n-32,  g=16  (data  is  2X  ovcrsampled) 

•  Vectorize  the  threshold  multiply  -  i.e.  Create  Vector  Y  =  1366  multiplies 

•  Vectorize  computation  of  Sl  and  of  Sr 

•  2  *  [  1 6  adds  to  get  started,  2  adds  (subtracts)  per  bin  after  that]  ~  4*  1 400  adds 

•  Zero  with  n+g  zeros  at  ends  for  startup  conditions 

•  Vector  compare  between  Y and  Sl.  and  Land  Sr 

•  Vectorize  the  logic  on  the  two  resulting  1-bit  w  ide  vectors 


Figure  26:  Standard  Cell-Average  Constant  False  Alarm  Rate  (CA-CFAR)  Detector 


5.5  Baseline  Radar  to  Processor  Interface 
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Figure  27  shows  a  data  flow  representing  a  baseline  radar  to  processor  interface.  Listed 
below  are  the  main  parameters 

•  ADC’s  run  at  3  Msps 

•  “Pulse”  is  4  ms  -  12k  record 

•  Detection  Mode 

•  1  Hz  PRF 

•  ADC  feeds  ARM  processor 

•  Imaging  Mode 

•  250  Hz  PRF  (=4  ms  PRI) 

•  Transmit,  receive,  and  processing  is  continuous 

•  ADC  feeds  FPGA-1 

Worst  case  analysis  list  is  captured  below: 

•  Target  is  at  the  maximum  range 

•  There  arc  a  series  of  targets  back-to-back  -  requires  continuous  imaging 

•  1 600  radar  records  -  6.4  seconds  per  image 

•  1 00  meter  range-swath  in  the  records 

•  250  Hz  PRF  (4ms  PRI) 

•  ARM  and  FPGA  have  memory  management  to  handle  continuous  imaging  on  multiple 
targets  (multiple  ranges  with  a  caravan  of  baek-to-baek  targets). 

12k  x  32  record  every  4ms 
while  imaging 


Figure  27:  Baseline  Scenario 


5.6  FPGA  processing 

To  have  zero-overhead  comer  turns  and  maximized  the  FPGA  logic  utilization,  there  arc  6 
memories,  each  with  its  own  bus,  mapped  to  the  pipelined  processing  stages.  Each  memory 
is  required  to  hold  the  entire  data  set.  In  order  to  do  comer  turns  at  the  maximum  rate  of  the 
pipelined  memory,  we  write  rows  and  read  columns  4  vectors  at  a  time.  We  describe  the 
processing  and  flow  of  data  through  the  6  memories  using  the  colored  chart  in  Figure  28. 
Eaeh  eolor  represents  a  pipelined  process  between  buffers  (i.e.  comer  turns).  The  goal  is  to 
maximize  FPGA  non-clockcd  and  power-down  durations.  FPGA-2  (Purple,  Brown,  and 
Blue  is  powered  down  while  FPGA-1  (Green)  is  processing. 
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Write-A 

Rows 

Ap-3 

Wnte-B 
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Ap-2 

Figure  28:  Data  flow  pipeline  through  the  6  memories 


5.6.1  FPGA-1  1-D  pulse  compression  processing 

Recall  how  dc-ramp  data  from  stretch  processing  is  generated,  as  illustrated  in  Figure  8.  The 
pulse  compression  processing  required  to  convert  the  data  into  the  equivalent  of  w  hat  would 
have  been  received  from  an  impulse  is  illustrated  in  Figure  29.  This  figure  shows 
background  subtraction,  filtering  and  decimation  of  the  ADC  samples,  the  Fourier  transform 
including  required  zero-padding  for  upsampling,  and  the  de-skew  operation.  Note  that  all 
these  operations  are  done  “at  the  same  time”,  that  is,  in  a  pipeline  that,  once  full,  provides 
output  data  in  real-time  with  the  input  data.  Figure  30  shows  operation  counts  for  the  various 
blocks.  The  Filter/Decimate  block  dominates  the  processing  load  and  sets  the  timeline. 
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Figure  29:  FPGA-1  Pulse  Compression  Flow 
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Figure  30:  FPGA-1  Pulse  Compression  Operation  Count 


IFwhcn  a  trajectory  estimator  is  added,  the  output  data  can  be  put  in  a  “spotlight”  mode, 
centered  on  the  target.  Figure  31  illustrates  spotlight  reformatting,  where  the  range  to  the 
first  point  in  a  record  is  adjusted  to  keep  the  target  in  the  center  of  a  short  record.  Basically, 
the  centroid  of  a  bright  point  is  maintained  in  the  center  of  the  collection  window,  and  the 
range  shift  performed  on  each  record  to  keep  it  centered  is  captured  so  that  later  stages  know 
how  to  interpret  the  data.  The  per  pulse  SNR  is  high  enough  to  do  single  pulse  detection  and 
to  allow  correlations  between  records  to  provide  the  time  lags  (and  phase  rotations)  required 
to  estimate  the  trajectory.  Trajectory  estimation  is  also  allows  computation  of  the  Stolt 
interpolator  and  phase  multiply  matrix. 

To  illustrate  the  data  reduction,  assume  we  want  an  18m  window  of  data  surrounding  the 
target.  Assume  a  range  bin  is  7.5  cm.  Then  the  number  of  output  range  bins  is  18m/7.5cm  = 
240  points  rather  than  the  1366  points  we  started  with. 

In  addition  to  allowing  computation  of  the  Stolt  interpolator  and  phase  multiply  matrix, 
trajectory  estimation  also  allows  presum  of  spotlight  centered  range  records.  This  pre¬ 
summing  would  drastically  cut  the  amount  of  data  that  would  need  to  be  exfiltrated  and/or 
processed  by  FPGA-2.  For  example,  the  aperture  for  a  target  at  maximum  range  contains 
1600  records.  If  we  can  presum  by  a  factor  of  8,  the  output  could  be  reduced  to  200  records. 
Order  3H  log2(H)  reduction  in  total  processing  from  could  be  gained  by  presuming  by  H.  In 
other  words,  the  computational  cost  of  performing  spotlight  reformatting  and  trajectory 
estimation  is  small  compared  to  the  savings  it  could  achieve.  Figure  32  depicts  spotlight 
processing  and  trajectory  estimation,  where  the  output  record  has  been  reduced  to  200  by 
240  points. 
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Figure  32:  FPGA-1  Pulse  Compression  Operation  Count 


5.6.2  FPGA-2  2-D  cok  Image  Former 

Figure  33  shows  the  major  processing  blocks  in  both  FPGAs  where  the  coloring  matches  the 
pipeline  sequencing  shown  in  Figure  28 — green  for  FPGA-1  and  purple,  brown,  and  blue  for 
FPGA-2.  Figure  34  is  a  timing  diagram  showing  how  the  various  functions  arc  done  in 
parallel  and  how  the  8k  FFTs  dominate  the  time  line.  Figure  35  shows  how  the  metric 
computations  arc  so  minor  they  can  be  ignored  for  the  purposes  of  power  calculations. 
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Figure  33:  Wavenumber  Migration  Processing  Partitioned  Around  Corner-Turns 
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Figure  34:  Example  of  Timing  on  FPGA-2  During  Image  Formation 
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Figure  35:  Metrics  Negligible  Processing  Load 


5.7  Power  Consumption  Analysis 

Summary  power  numbers  for  FPGA-1  arc  shown  in  Figure  36.  Summary  power  numbers  for 
FPGA-2  arc  shown  in  Figure  37.  The  numbers  are  based  on  back-to-back  targets 
continuously  appearing  at  the  farthest  range.  At  closer  ranges  the  power  will  be  lower 
because  an  aperture  has  fewer  pulses.  While  the  FPGA-2  static  power  is  drastically  reduced 
by  power  cycling,  the  remaining  static  power  dissipation  is  caused  by  the  configuration  time 
each  time  the  FPGA  is  powered  up.  FPGA-1  static  power  is  high  because  there  is  not 
enough  time  to  power  it  down  given  the  time  it  takes  for  configuration  before  it  must  he 
ready  to  operate  again.  Using  a  process  like  Altcra’s  “Hardcopy™”  would  eliminate  the 
configuration  time  and  thereby  eliminate  this  configuration  time  power,  and  in  the  case  of 
FPGA- 1 ,  the  time  it  remains  powered  up,  when  it  could  be  powered  down. 
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Figure  36:  FPGA-1  Power  Summary  Calculation 
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igure  37:  FPGA-2  Power  Summary  Calculation 


Another  option  to  reduce  power  further  would  be  to  add  a  smaller  FPGA  to  just  manage  a 
buffer  between  the  radar  and  the  A/B  RAM.  The  task  would  be  storing  N=64  records 
(64*4ms=256ms)  and  bursting  them  to  RAM-A/B.  FPGA-1  could  now  operate  in  a  burst 
mode  as  follows: 

•  60ms  wakc/config,  10ms  to  process  64  pulses  &  go/stop  test,  Power  off  for  186  ms 

•  70/256  duty  cycle  would  reduce  static  power  by  1/4  --  i.c.  from  150mW  to  41  mW 
Summary  power  numbers  for  the  memory  is  shown  in  Figure  38.  Summary  power  numbers 
for  the  ARM  processor  is  shown  in  Figure  39. 
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Figure  38:  Memory  Power  Summary  Calculation 
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Figure  39:  ARM  Power  Summary  Calculation 


5.8  Digital  Hardware  Power  Summary 

Figure  40  sums  the  power  drawn  by  the  individual  parts. 
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ation 


We  have  shown  that  the  processing  can  be  done  in  a  small  ARM  processor,  2  FPGAs,  and 
12  cell-phone  RAM  chips.  We  have  shown  that  the  power  drawn  in  the  detection  mode  is  ~5 
mW.  For  imaging,  we  have  shown  that  a  total  power  draw  of  <200  mW  in  a  “worst  case” 
scenario  where  a  maximum-range  target  is  in  scene  continuously.  For  ~$200k  NRE  (non 
recurring  engineering  cost),  the  “HardCopy™”  process  can  be  applied  to  harden  FPGA-1, 
resulting  in  a  total  power  draw  of  <100mW. 
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6  Antenna  Design 

Wide  bandwidth  antennas  are  often  dispersive.  Figure  41  illustrates  most  of  the  classie 
wideband  antenna  designs.  The  fat  dipole  and  the  TEM  horn  are  the  only  ones  on  this  chart 
that  are  not  dispersive.  In  order  to  obtain  high  resolution  SAR  and  ISAR  imagery,  the 
antenna  must  have  low  dispersion,  and  have  low  dispersion  across  the  wide  integration  angle 
of  the  SAR  or  ISAR  system.  If  the  antenna  does  not  inherently  have  low  dispersion,  then  the 
dispersion  ean  be  theoretically  be  processed  out  by  deconvolution — but  only  if  the 
dispersion  is  consistent  over  the  wide  integration  angle  of  the  SAR  or  ISAR.  Figure  42 
illustrates  this  deconvolution  processing.  In  this  ease,  the  antenna  has  a  transfer  function  of 
H(w)  and  the  deconvolution  filter  Q(w)  =  1/H(w).  The  key  is  that  the  resulting  impulse  must 
be  stationary  relative  to  the  antenna  pointing  angle.  In  other  words,  the  two  impulses 
pictured  on  the  oscilloscope  must  not  move  when  the  antenna  is  rotated  about  its  phase- 
center. 
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Figure  41:  Wide  bandwidth  antennas  are  often  dispersive 
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Figure  42:  Setup  to  measure  angle  dependent  dispersion 

Since  the  FCC  adopted  its  UWB  regulations,  significant  effort  has  gone  into  building  small 
non-dispersive  ultra  wide  bandwidth  antennas.  Figure  43  is  a  picture  of  a  set  of  these.  The 
third  from  the  right  on  the  top  row  is  a  mcanderline  antenna.  The  majority  arc  variations  of  a 
fat  dipole  or  “bowtic”  antenna  with  different  feeds.  Some  are  electrically  small  dielectric 
loaded  elements  that  sit  on  a  large  “ground  plane.”  In  reality,  the  so-called  ground-plane  is 
actually  the  main  radiator  and  essentially  the  other  half  of  a  dipole  formed  with  the 
dielectrically  loaded  clement.  Indeed,  a  major  problem  with  these  antennas  was  the  cable 
radiation  caused  by  the  antenna  using  “ground”  currents  on  the  shield  of  the  coax.  A  pair  of 
dielectrically  loaded  elements  connected  with  a  balun  resulted  in  a  small  radiator  with  good 
efficiency  and  low  cable  radiation.  The  two  antennas  on  the  lower  left  can  be  thought  of  as  a 
folded  slot,  or  a  top  loaded  monopole,  again  with  different  feed  structures.  These  antennas 
had  the  highest  and  flattest  gain  of  all  the  antennas  pictured,  as  well  as  reasonably  low  shield 
currents  when  the  coaxial  feed  connection  was  placed  on  the  side,  rather  than  bottom  of  the 
antenna  (as  shown  on  the  bottom  left  antenna). 
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Figure  43:  UWB  antennas  developed  since  the  FCC  issued  UWB  rules 

Figure  44  is  an  expanded  view  of  this  antenna  showing  several  feeds.  This  antenna  has  the 
following  features  that  are  important  for  this  application: 

•  It  is  electrically  small  ~  14  X  at  lowest  frequency,  where  the  return  loss  is  only  6  dB 

•  It  is  Ultra  Wide  Bandwidth  -  The  -6  dB  bandwidth  is  about  the  same  as  the  center 
frequency.  Figure  45  shows  the  transfer  function  of  the  antenna. 

•  It  can  have  a  tuning  clement  that  creates  a  narrow  notch  in  the  spectrum  if  needed  to 
reject  interference  (US  Patent  6,590,545  and  6,914,573  and  7,352,333) 

•  It  has  a  null  in  the  plane  of  the  antennas  that  will  enhance  isolation  between  Tx  and  Rx 
antennas.  This  is  important  because  Tx  Rx  coupling  is  effectively  a  “jamming”  term  that 
stresses  the  dynamic  range  of  the  radar. 

•  It  can  be  conformed  to  a  curved  shape  if  necessary  —  it  docs  not  have  to  be  planer. 

•  It  has  a  broad  beam,  which  is  needed  to  support  wide  SAR  and  1SAR  integration  angles. 
Its  peak  gain  is  only  about  5  dB  on  boresite  at  the  peak  frequency. 

•  It  has  low  dispersion  -its  impulse  response  has  low  time  sidclobcs.  This  is  important  in 
that  battery  draining  processing  is  not  required  to  deconvolve  the  antenna  response. 
Down-range-resolution  will  be  close  to  ideal  without  deconvolution.  If  deconvolution 
corrections  are  used,  they  will  be  mild  and  arc  not  likely  to  cause  serious  artifacts. 

Figure  46  shows  the  impulse  response  of  a  pair  of  antennas  -  in  other  words,  the  round- 
trip  radar  response.  This  data  is  the  derivative  of  step  response  taken  with  a  TDR 

•  The  impulse  response  is  stable  across  wide  integration  angles.  This  is  important  because 
the  1SAR  integration  gain  and  eross-range-resolution  will  be  close  to  ideal.  Figure  47 
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shows  the  impulse  from  several  orientations  overlaid  to  illustrate  the  stability.  This  data 
is  the  real  part  of  the  inverse  Fourier  transform  (IFFT)  of  S21  measured  with  a  vector 
network  analyzer  (VNA). 


•  It  can  be  backed  by  a  reflector.  The  first  test  antenna  used  in  the  MTR1  field  tests  used  a 
transmit  and  reecive  element  baeked  with  an  aluminum  plate.  The  baek  plate  enhanced 
the  midband  gain  and  gave  the  antenna  a  good  front  to  back  ratio. 

•  A  cavity  backed  version,  with  a  modified  shape,  is  also  viable.  This  version  was  used  in 
MTRl’s  second  field  test.  The  cavity  backed  version  improved  the  Tx-Rx  isolation. 


Max  radiation  normal  to  plane 

Figure  44:  Folded  Slot 


Figure  45:  Facc-to-Facc  Two-Antenna  S21  (nose-to-nose  S21,  NOT  gain) 
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Figure  47:  Tight  Impulse  Response  Versus  Orientation 

To  further  understand  the  antenna  characteristics,  chamber  measurements  were  performed 
on  a  scaled  antenna  -  scaled  up  in  frequency  by  a  factor  of  5.7.  Figure  48  shows  the 
orientation  nomenclature  used  in  plotting  the  results.  Figure  49  shows  a  co-pol  E-planc  and 
H-planc  frequency  vs  angle  energy  plot.  Note  the  natural  isolation  in  the  plane  of  the 
antenna  at  90  and  270  degrees  in  the  H-planc  cut  data.  Note  also  the  broad  flat  radiation 
pattern  in  the  +/-  30°  from  borcsite  in  both  elevation  &  in  azimuth.  Figure  50  shows  a  cross- 
pol  E-planc  and  H-planc  frequency  vs  angle  energy  plot.  Figure  51  shows  beam  patterns  at 
4,  5,  and  8  GHz  which  scale  to  700MHz,  877MHz,  and  1.4  GHz.  These  patterns  show 
excellent  wide  angle  wide  bandwidth  imaging  capability. 
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Figure  48:  Measurement  Orientation  Nomenclature 
H-Plane  Cut  E-Plane  Cut 


Figure  49:  Co-pol  E-planc  and  H-plane  Frequency  vs  Angle  energy  plots 


TR-1 10901  Applied  Signals  Intelligence,  Inc.  FOR  OFFICIAL  USE  ONLY  (FOUO)  Page  45  of  49 


Subagreement  #100541Z1 


Riverine  Imagin 


September  30,  2011 


H-Plane  Cut 


E-Plane  Cut 


Figure  51:  Beam  Patterns  Sealed  to  700  MHz,  877  MHz,  &  1.4  GHz 
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Two  antennas  were  fabricated  for  MTRI.  Figure  52  is  a  picture  showing  the  first  antenna 
being  mounted  to  a  tree.  It  consisted  of  a  pair  of  the  folded-slot  antennas,  one  for  transmit 
and  the  other  for  receive,  mounted  in  front  of  a  reflector  plate.  Figure  53  shows  a  picture  of 
the  second  antenna  with  the  cavity  backing.  The  cavity  backing  that  improves  the  isolation 
between  the  Tx  and  Rx  antenna  elements.  The  antenna  was  slightly  larger  than  the  first  to 
move  the  bandwidth  to  lower  frequencies  that  would  penetrate  wet  wood  better. 


Figure  52:  First  Antenna  Delivered  To  MTRI 
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Figure  53:  Second  Antenna  Showing  Cavity  Backing 

Figure  54  shows  data  collected  from  the  cavity  backed  antenna  delivered  to  MTR1  for  its 
second  field  test.  The  reflection  coefficient  measured  with  a  TDR  shows  a  relatively  smooth 
taper  from  50  ohms  to  an  open  that  approximates  a  Klopfenstein  taper.  As  a  result,  the  plots 
of  return  loss  and  VSWR  show  the  antenna  match  improves  with  frequency  and  remains 
good  as  the  frequency  goes  up.  The  nose-to-nosc  impulse  response  shows  low  sidelobes,  and 
the  transfer  function  plot  shows  +/-  1  dB  of  ripple  through  the  passband  of  the  antenna.  This 
ripple,  and  the  sidelobes  in  the  impulse  response  could  be  improved  by  loading  the  cavity, 
but  at  the  expense  of  reduced  gain. 

In  conclusion,  an  antenna  small  enough  to  fit  the  intended  application  was  fabricated  and 
proved  have  all  the  performance  metrics  needed  to  perform  the  high  resolution  UWB  ISAR 
function. 
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reflection  coeff 


Sit  Return  Loss  dB 


VSWR 


2-antenna  round-trip  impulse  response 


2-antenna  round-trip  transfer  function 


Figure  54:  Trimmed  ASI390-T2.5-TG1 .8 
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